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Summary 


Love waves and Rayleigh waves from eight earthquakes recorded 
at Hallett Station, Scott Base and Mirny have been analysed and the 
dispersion compared with that predicted by theory for simple model 
crusts. The average thickness of the crust in eastern Antarctica is found 
to be about 35 km, as is typical of continents, whereas Marie Byrd Land 
with an average thickness of about 25 km cannot be regarded as truly 
continental. Love wave dispersion indicates that the thickness of the 
solid crust in the oceanic regions surrounding Antarctica varies from 
about 5km to 1okm, the smaller values being associated with the 
deeper basins. It is shown that the determination of oceanic crustal 
thickness from Rayleigh wave dispersion is in general subject to large 
uncertainties, nor can one usually rely on values of the thickness of 
unconsolidated bottom sediments obtained by this means. 


1. Introduction 


The broad structure of the land beneath the Antarctic ice cap can be determined 
only by indirect means. The ice cap is of continental size but this does not imply 
a single co-extensive land mass. Geological evidence points to the existence of a 
continental shield in eastern Antarctica whilst western Antarctica appears to be 
characterized by young folded belts, but the extent of these two geological pro- 
vinces is uncertain. Much of the coastline is devoid of exposed rock and little is 
known of the region between the Weddell Sea and the Ross Sea, separating eastern 
from western Antarctica. Seismic reflection work in Antarctica has revealed that 
the rock surface lies below sea-level at many places in the interior, and often at a 
sufficient depth to remain below sea-level if the ice were to melt, even allowing for 
consequent isostatic uplift. 

The installation of seismograph stations in Antarctica for the International 
Geophysical Year has provided an opportunity for determining the main features 
of the Antarctic crust by a study of guided earthquake waves. In a preliminary 

? 289 





Editors 
A. H. Cook, M.A., Ph.D., F.RAS., F.G.S. —_T. F. Gaskell, 


: M.A., Ph.D., F.RA-S. 
Natisnal Phoceal spite <. British Petroleum Company, 


Teddington. 


The GonplyslalSeornal in setpil See ts gubiiaadad oso gegliyeten 
and related subjects. ee ee 
Editors are accepted for publication. 


The aims of the Editors are to papers of a high standard, to give them 
« ide innencluindl diaeloiny cade mate aed eaeny publication. 


Book reviews and reports of a a2 re 
ee oe ON et ee ee Oa ee 


Notice to Authors 
Authors of papers offered for publication need not be Fellows of the Society. 


ribet aiacauiliod cs waite 25 reais Gaus saben ac 
Contributors are entitled to receive 25 reprints numbers are 
obtainable at moderate rates, ¢,g. 300 copies of a 12-pgs reprin rint oe. plates) 
cost £11 14s. (US, $35). Spec Special arrangements can ai goer 
requiring reprints contributions with covers, special pagination etc., for 
issue in their own series of publications. 


Attention is drawn to the recommendations to authors printed on the inside 
of the back cover, 





Volume 3 No. 4 December 1960 


The following contributions have been: accepted for publication 
in the above issue: 


R. J. Donato, Experimental investigation on the properties of Stoneley 


waves. 
T. Haruerton, ‘Mierceciatia at Scott Base. 


B. A. Bo.t, The revision of earthquake epicentrés, focal depths and origin- 
times, using a high-speed computer. | 


E. InvinG, Palaeomagnetic directions and pole positions, Part IT. 


R, W. Wiitmore & A. M. BANCROFT, The time-term —— to refraction 
seismology, 


M. J. P. Mosca, Refleson and refraction of plane ate waves at «plan 
boundary between aeolotropic media. 


J. C. Jarcer & R. F. Toran, Geapjeics in Aaya: 





The Geophysical - Journal 
of the 


ROYAL ASTRONOMICAL SOCIETY 


Vol. 3 No. 3 September 1960 





Thickness of the Earth’s Crust in Antarctica and the 
Surrounding Oceans 


F. F. Evison, C. E. Ingham, R. H. Orr and J. H. Le Fort 


(Received 1960 February 10) 


Summary 


Love waves and Rayleigh waves from eight earthquakes recorded 
at Hallett Station, Scott Base and Mirny have been analysed and the 
dispersion compared with that predicted by theory for simple model 
crusts. The average thickness of the crust in eastern Antarctica is found 
to be about 35 km, as is typical of continents, whereas Marie Byrd Land 
with an average thickness of about 25 km cannot be regarded as truly 
continental. Love wave dispersion indicates that the thickness of the 
solid crust in the oceanic regions surrounding Antarctica varies from 
about 5km to 1okm, the smaller values being associated with the 
deeper basins. It is shown that the determination of oceanic crustal 
thickness from Rayleigh wave dispersion is in general subject to large 
uncertainties, nor can one usually rely on values of the thickness of 
unconsolidated bottom sediments obtained by this means. 


1. Introduction 


The broad structure of the land beneath the Antarctic ice cap can be determined 
only by indirect means. The ice cap is of continental size but this does not imply 
a single co-extensive land mass. Geological evidence points to the existence of a 
continental shield in eastern Antarctica whilst western Antarctica appears to be 
characterized by young folded belts, but the extent of these two geological pro- 
vinces is uncertain. Much of the coastline is devoid of exposed rock and little is 
known of the region between the Weddell Sea and the Ross Sea, separating eastern 
from western Antarctica. Seismic reflection work in Antarctica has revealed that 
the rock surface lies below sea-level at many places in the interior, and often at a 
sufficient depth to remain below sea-level if the ice were to melt, even allowing for 
consequent isostatic uplift. 

The installation of seismograph stations in Antarctica for the International 
Geophysical Year has provided an opportunity for determining the main features 
of the Antarctic crust by a study of guided earthquake waves. In a preliminary 
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notice of the work described in this paper it was shown that the eastern Antarctic 
crust in the region from Victoria Land to Wilkes Land is about 35 km thick 
(Evison, Ingham & Orr 1959). This finding provided definite geophysical support 
for the presumption of a land continent in Antarctica, since the continents are 
characterized by a crustal thickness of about 35 km in contrast to thicknesses of 
5-15 km in oceanic regions. Similar studies of the dispersion of Love and Rayleigh 
waves have also been carried out by Press & Dewart (1959), from seismograms 
obtained at Wilkes Station, and whilst they confirm the continental thickness of the 
area mentioned above, they report a much thinner average for paths through other 
sectors of Antarctica. This latter result is partly explained by the present study. 

Hallett Station and Scott Base are well placed for distinguishing between 
eastern and western Antarctica, and separate determinations have been obtained 
for the crustal thickness in these two major sectors. Records from Mirny have 
been used mainly for measuring the dispersion of waves travelling across the 
South Indian Ocean. It appears that for most of the paths along which Press 
and Dewart found the average thickness to be less than continental much of the 
deficiency occurs in western Antarctica. The present results indicate, however, 
that the crust is continental throughout eastern Antarctica. 
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Fic. 1.—Wave paths, showing thickness of crust inferred for each segment. 
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2. Observations 


Details of the eight earthquakes used in this study are listed in Table 1. 
Epicentres and origin times were taken from the monthly bulletins prepared by 
BCIS and USCGS. Earthquakes were chosen for which the wave paths would 
give good coverage of Antarctica, as illustrated in Figure 1. Epicentres close to 
Antarctica were preferred so that the contribution of the Antarctic segment of the 
path to the total dispersion would be as large as possible. Details of the seismo- 
graphs that provided the records used in this study are as follows: 

Hallett Station: Columbia seismograph, 3 components 
To = 158, Ty = 758. 

Scott Base: Benioff seismograph, vertical and E-W, 
To = 0-68, Ty = 258. 

Mirnv’ Kirnos seismograph, 3 components, 
To = 12°58, Ty = 1°28. 

The preference for epicentres close to Antarctica results in shorter paths being 
used than is usually the case in studies of this kind. Consequently the degree of 
dispersion on some records was rather limited. A majority of the records studied, 
however, came from the long period Columbia seismograph at Hallett Station, 
and a few of these are reproduced in Figure 2. 
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Fic. 2.—Rayleigh waves recorded at Hallett Station for earthquakes of 
1957 Sept. 4 (a and b), 1957 Oct. 2 (c), and 1957 Oct. 20 (d). 
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3- Method of analysis 


The virtual absence of seismic activity in Antarctica makes it necessary to rely 
on waves that have first passed through some neighbouring oceanic region, and 
perhaps through some other continent as well. The recorded dispersion is there- 
fore composed of two or more separate effects, of which one is of primary interest 
to us. The usual method of dealing with such composite wave paths is to assume 
that the dispersion due to an extraneous segment conforms with an empirical 
curve that has been derived for a known region of normal continent or ocean. 
When a correction has been made for this segment the remaining dispersion may 
also be assessed by comparison with the appropriate empirical curve. This is the 
method used by Press & Dewart (1959) in their analysis of Love and Rayleigh 
waves across Antarctica. The advantage of the method is that one does not have 
to postulate a simplified model crust and assume values for the various seismic 
parameters. 

In the present study a different method has been adopted, in which the observa- 
tions are compared with theoretical rather than empirical results. Theoretical 
dispersion curves have been computed from equations given by Stoneley (1948) 
for Love waves, by Lee (1932) for continental Rayleigh waves and by Jardetzky & 
Press (1953) for oceanic Rayleigh waves. The advantage of this method is that any 
departure from normal crustal thickness is readily accommodated. One finds that 
where a well-recorded dispersive wave has traversed a composite path the analysis 
reveals a value of crustal thickness for each segment, and surprisingly often there 
is a unique combination of values that will give reasonable agreement between 
observation and theory. Thus in the present study the results for Antarctica do 


not depend on assumed values of crustal thickness for the adjacent regions. 

The models adopted for the purpose of computing theoretical dispersion curves 
consist simply of horizontal layers of uniform thickness, overlying a semi-infinite 
mantle. For the continental case a single crustal layer is postulated, and for the 
oceanic case a water layer and a single rock layer. The following seismic constants 


are assumed, incorporating the parameters used in Stoneley’s (1948) computations 
for Love waves: 


Water Crust Mantle 
Compressional velocity (km/s) 1°52 6-01 7°79 
Shear velocity (km/s) 3°47 4°50 
Crust/water Mantle/crust 
Density ratio 2°65 1°283 


Figure 3 shows the family of theoretical Love wave curves for selected thick- 
nesses of crust. Since Love waves are not affected by a fluid layer the dispersion 
for any path, whether continental or oceanic, is given by an appropriate curve of 
this family. Figure 4 shows dispersion curves for continental Rayleigh waves. 

In computing the theoretical dispersion of oceanic Rayleigh waves one has to 
assign a further parameter, namely the ratio of the thickness of the rock layer 
(Hj) to that of the fluid layer (Ho). ‘Table 2 shows values of U and T/Ho, where 
U is the group velocity and 7 the period, computed for H;/Ho = 1, 2 and 3. 
Dispersion curves for selected values of Hj are given in Figure 5. 

In a wave path that is made up of a number of segments let my be the length 
of the mth segment, expressed as a fraction of the whole path length, and let Up, 
be the group velocity of a particular period traversing that segment. Then the 
group velocity of the period taken over the whole path is U = 1/2p(mp/Un). 
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Thus we may construct theoretical dispersion curves for a composite path by assign- 
ing to each segment some curve from Figure 3, 4 or 5, making due allowance in 


Table 2 
Dispersion of oceanic Rayleigh waves 


M/Ho = 1 Hi/Ho = 2 ;/Ho = 3 

U T/Ho U T/Ho U T/Ho 
(km/s) (s/km) (km/s) (s/km) (km/s) (s/km) 
1°157 2°595 1°350 2°719 1°584 2°917 
1684 2°723 1°813 2°880 2°123 3°165 
2°139 2°856 2°421 3°065 2°642 3°445 
3°01! 3°138 3°183 3°474 3°244 4021 
3°689 4°195 3°747 5°179 3°776 6-309 
3°859 5°239 3°g900 6-986 3°905 8-831 
4°002 8-210 4°00! 11-908 3°998 15°442 


each case for the appropriate fractional length m,. In measuring these lengths we 
have adopted the usual convention that the boundary between land and ocean 
is defined by the 1000 fathom depth contour. In Figure 1 each wave path is 
shown divided into its segments; the captions to Figures 6-13 give the length of 
each segment as a percentage of the whole path. For simplicity, in a few cases 
a short stretch of water shallower than 1000 fathoms has been included in an 
oceanic segment. 

The recorded dispersion is represented by means of points on a graph; this is 
achieved by the usual procedure, in which a graph of arrival time against serial 
number is prepared for a suitable succession of recorded phases. The interpreta- 
tion of the observed dispersion results for a particular path is then a matter of find- 
ing those theoretical curves which, when combined in the manner described above, 
provide the best fit to the plotted data. 

The average depth of water along each oceanic segment has been calculated 
from bathymetric charts. This depth is useful, as will be seen later, in helping to 
resolve the indeterminacy that usually arises in interpreting the dispersion of oceanic 
Rayleigh waves. 


4- Choice of parameters 


The type of crustal model that one must be content to adopt in this kind of 
study, if the computations are to be performed in a reasonable time by desk calcu- 
lator, is undoubtedly very much simpler than any real crust. Consequently it is 
not to be expected that the seismic parameters appropriate to the model will cor- 
respond precisely with particular values that might be determined by other means, 
for example by a study of body waves from earthquakes or artificial explosions. 
Such values usually pertain to a limited region of the medium—thus refracted 
waves are apt to indicate the velocity near the upper surface of the refracting 
layer—and they cannot well be adopted for theoretical purposes until more complex 
models can be readily analysed. 

Even the simple models adopted here involve seven seismic parameters, apart 
from the ratio of fluid depth to crustal thickness that is required in the case of 
oceanic Rayleigh waves. Computations can be found in the literature for a variety 
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of parameter values, but a thorough study of the effect of each separate parameter 
upon the dispersion curve must await a systematic coverage of the problem by 
electronic computer. In the absence of any compelling reason for a more elaborate 
choice of parameters it has been thought desirable to use the same values for the 
continental and oceanic cases, and also for both Love and Rayleigh waves. 

The values of compressional velocity (6-01 km/s) and shear velocity (3-47 km/s) 
chosen for the crust are close to those preferred by Press & Ewing (1955) for the 
continental case. For the oceanic crust these authors adopt much greater values, 
especially in their Rayleigh wave computations. According to Hill (1957), how- 
ever, the evidence of seismic refraction work is that in water depths greater than 
4km the solid crust comprises on an average a layer 1°75 km thick with compres- 
sional velocity 4-6km/s and a layer 4-7km thick with a velocity 6-71 km/s. Our 
value of 6-01 km/s may be a little too low, but is not inconsistent with Hill’s pic- 
ture, even though the water depth for some paths is slightly less than 4 km. 

The compressional velocity 7°79 km/s chosen for the mantle is at about the 
lower limit of the range of values that have been measured in various parts of the 
world by means of refracted waves. This low value, together with the correspond- 
ing shear wave velocity (4-50km/s) gives a better fit in dispersion studies than do 
values such as 8-1km/s and 4°7km/s, which are more commonly indicated by 
refraction. This is a case where it seems inappropriate to attach more typical 
measured values of the parameters to our simplified models; conversely, it is not 
to be thought that a refraction survey would necessarily indicat- the same values 
as those adopted here. 

The general validity of the present choice of parameters may be tested by com- 
paring the resulting theoretical dispersion curves with empirical curves that have 
been obtained for paths of known type. Such empirical curves, as given by Press 
& Dewart (1959) for continental paths, are included in Figures 3 and 4; it will be 
seen that they agree quite closely with the theoretical curves for a crustal thickness 
of 35km, which is usually regarded as the average thickness for continents. 


5. Results 


Figures 6-13 show the plots of observed dispersion for each wave path and 
type of wave, together with the theoretical composite curves that appeared to give 
the best fit. Two or three such composite curves are presented in each figure so 
as to give some impression of the sensitivity of the method. In most cases the fit 
that would be obtained by any other interpretation is distinctly inferior, as may 
readily be checked by reference to the theoretical curves in Figures, 3, 4 and 5, 
taking into account the percentage of path contributed by each particular 
segment. 

The most direct evidence that the thickness of the crust in eastern Antarctica 
averages about 35 km is given by the Love wave analyses in Figures 6 and 7. In 
both cases a wide range of periods was recorded and the analyses yielded inde- 
pendent determinations of thickness for the oceanic segments. The result was 
confirmed by analyses of Rayleigh waves from the same two earthquakes, as 
shown in Figures 8 and g. Further confirmation from another earthquake is 
given in Figure 10a. A valuable determination of the same thickness is that of 
Figure 11a, since the path crosses Antarctica diametrically and the oceanic segment 
is small. Results for Rayleigh and Love waves recorded at Mirny for an epicentre 
in south Bolivia are included (Figures 10b and 12c) as offering an interpretation 
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consistent with that already put forward, although only a small range of periods 
was recorded. 

A much thinner average crust is indicated for western Antarctica. Four wave 
paths have been studied, all traversing Marie Byrd Land, and although the 
observations did not provide so definite a determination as in the case of eastern 
Antarctica, it appears that the average thickness in this region is about 25 km. 
The simplest demonstration of this result is given by Figure 12b, which gives the 
dispersion of a Love wave from an epicentre in Drake Passage, recorded at Hallett 
Station after traversing some 1 500 km of ocean and 3 500km in western Antarctica. 
A glance at the family of theoretical Love wave curves (Figure 3) will serve to show 
that these observations do not allow of a crustal thickness much greater than 
25km for the Antarctic segment. 

The same conclusion arises from the Rayleigh wave analysis of Figure 11b for 
an epicentre in south Bolivia. This path traversed peninsular South America for 
43 per cent of its length, and by allowing a crustal thickness of 30km for this 
segment a very satisfactory fit has been obtained. This value seems rather thin 
for such a region and needs further verification; but if a thicker crust were assumed 
for the South American segment the western Antarctic crust would have to be 
taken thinner, and then also the degree of fit in Figure 11b would be inferior. 
Figure 13a also indicates a thickness of about 25km for western Antarctica 
if we assume the normal continental value of 35km for the Brazilian shield. 
Further Rayleigh wave data of a more limited kind are presented in Figures 
13b and 11c, where they are shown to give a consistent solution for western 
Antarctica. 

The thickness of solid crust beneath the surrounding oceans can be uniquely 
determined only from Love waves, which are unaffected by the water layer; 
furthermore, an accurate determination depends on the shorter periods. As we 
have seen, Figure 6 indicates a thickness of no more than 5 km for the southwest 
Indian Ocean, and Figure 7 indicates about 10okm for the southeast Indian 
Ocean. Values close to 10okm have been obtained from Love waves for Drake 
Passage (Figure 12b) and for the region northeast of the Ross Sea (Figure 
12a). The results shown in Figure 12a also involve a short segment of path across 
the eastern Ross Sea, which appears to have a crust of almost continental thick- 
ness, and thus to belong to eastern rather than western Antarctica. 

These values of crustal thickness for oceanic regions, obtained from record- 
ings of Love wave dispersion, appear to be related to the depth of the ocean much 
in the manner that one would expect on the principle of isostatic balance. As 
Table 1 shows, the water was 4.50km deep on an average for the path where a 
5km crust was indicated (Figure 6) and about 3-75 km deep for the paths with 
a crust of about 10km thick (Figures 7, 12a, 12b). This relation has been used in 
the interpretation of oceanic Rayleigh waves, which otherwise, as will be dis- 
cussed further below, do not yield an unambiguous solution. The major uncer- 
tainty has been resolved by taking Hi/Ho = 2 for the water depths from 3-71 km 
to 3-98km, and H,/Ho = 1 for depth of 4-18km to 4-50km (see Table 1). Thus, 
although a close fit has been achieved in most cases between the theoretical and 
observed Rayleigh wave dispersion for oceanic segments, this part of the analysis 
is essentially dependent upon the Love wave results and the estimates of water 
depth. The method is capable of further refinement by using fractional values of 


H,/Ho, but this is not warranted by the limited quantity of data used in the present 
study. 
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Fic. 10.—Rayleigh wave dispersion: 
(a) Indian Ocean to Hallett (60 per cent Indian Ocean at os 40 per 
cent E. Antarctica at 30-40km). Components: x E-W 
(b) S. Bolivia to Mirny (27 per cent S. America at 30km; 35 et cent 
S. Atlantic at 9 km; 38 per cent E. Antarctica at 30-40 km). Components: 
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Fic. 11.—Rayleigh wave dispersion: 
(a) Bouvet Island to Hallett (28 per cent S. Atlantic at 5km; 72 per 
cent E. Antarctica at 30-40km). Component: + Z. 
(b) S. Bolivia to Hallett (43 per cent S. America at 30km; 26 per cent 
S. Pacific at 5-5 km; 31 per cent. W Antarctica at 20-25 km). 
Component: x Z. 
(c) Drake Passage to Hallett (30 per cent Drake Passage at 8km; 70 per 
cent W. Antarctica at 20-25km.) Components: eE-W; © Z. 
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Fic. 12.—Love wave dispersion: 
(a) S. Pacific to Scott (80 per cent S. Pacific at 10-124km; 20 per cent 
Ross Sea at 30km). Component: e E-W 
(b) Drake Passage to Hallett (30 per cent Drake Passage at 9km; 70 per 
cent W. Antarctica at 20-25 km). Component: + E-W 
(c) S. Bolivia to Mirny (27 per cent S. America at 30km; 35 per cent 
S. Atlantic at 9km; 38 per cent E. Antarctica at 30-40km). Compo- 
nent: © N-S. 
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Fic. 13.—Rayleigh wave dispersion: 

(a) Mid-Atlantic to Hallett (9 per cent Atlantic at 9km; 28 per cent 
S. America at 35km; 23 per cent S. Atlantic at 9km; 40 per cent W. 
Antarctica at 20-25km). Components: e N-S; © Z. 

(b) S. Bolivia to Scott (47 per cent S. America at 30km; 19 per cent 
S. Pacific at 5km; 34 per cent W. Antarctica at 20-25km). Compo- 
nent: + Z. 
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Several wave paths in Figure 1 make rather small angles with the 1000 fathom 
depth contour. In such cases theory shows that the actual path is deflected from 
the great circle because of refraction at the boundary. The case where this effect 
would be expected to be most marked is the path from south Bolivia to Hallett 
Station. A careful study of this path has shown, however, that the result arrived 
at in Figure 11b would not be substantially altered by allowing for refraction. 


6. Indeterminacy of oceanic Rayleigh wave solutions 


The approximate thicknesses of oceanic crust shown in Figure 1 are com- 
patible with the observed dispersion of Rayleigh as well as Love waves, but could 
not have been ascertained from Rayleigh waves alone, except between rather wide 
limits. It has been explained that the Love wave results support the idea of an 
inverse relation between water depth and crustal thickness, as is required by the 
principle of isostasy. The Rayleigh wave solutions have been obtained with the 
help of this added condition. In the Rayleigh wave theory, if the thickness of the 
fluid layer could be taken as equal to the known depth of water, an analysis of 
Rayleigh wave dispersion would yield a unique value of crustal thickness; the 
analysis would then consist of finding the value of H/Ho that gave the best fit. 
It is generally found, however, that better agreement with other lines of evidence 
can be achieved by assuming a somewhat thicker fluid layer, the additional thick- 
ness being ascribed to unconsolidated bottom sediments. 

Numerous studies appear in the literature in which estimates of the thickness 
of the oceanic crust have been based upon Rayleigh wave dispersion. Useful sum- 
maries may be found in review articles by Press & Ewing (1955) and Hill (1957). 
The claim that a thickness of about 5 km is normal for all oceanic areas is largely 
based on Rayleigh wave data. But these studies do not seem to take account of the 
fact that very similar dispersion may be produced by widely different crustal 
thicknesses, as readily follows from the theory if we allow that variations may 
occur in the ratio of crustal thickness to fluid depth (H)/Hp). 

The point is illustrated in Figure 14, which gives three different interpretations 
for a Rayleigh wave recorded at Mirny from the earthquake of 1957 August 4, the 
path being virtually oceanic throughout (Figure 1). The three theoretical curves 
are computed for the same seismic parameters, but H)/Hp is taken as 1, 2 and 3 in 
turn. To the accuracy that can reasonably be claimed for dispersion studies there 
is little to choose among the three solutions, yet the crustal thicknesses indicated 
are respectively 5} km, 9$km and 12km. (The Love wave analysis suggests that 
the first of these solutions is nearest to the correct one.) Naturally any other value 
within this range could equally well be deduced by taking the appropriate fractional 
ratio of the thicknesses. The range of possible solutions might be less in cases 
where the observations extend over a greater range of velocity; but the later por- 
tions of most records bear signs of interference, indicating that energy has been 
propagated along paths other than the direct one, and hence making the analysis 
less reliable. 

Rayleigh waves have also often been invoked as a means of estimating the thick- 
ness of the unconsolidated bottom sediments. Inferred values are usually in the 
range 0-4-1'2 km; thus the effect of the sediments is a second-order one. But the 
actual value that is obtained in any particular case is wholly dependent on the as- 
sumed ratio of crustal thickness to fluid depth. Thus in Figure 14 the three inter- 
pretations give the thickness of the fluid layer as 5} km, 43km or 4km. The 
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average water depth for this path as determined from the bathymetric map is 
4'18km. Hence one inferred value of the sediment thickness is 1-07 km, another is 
o*57km, while the third interpretation would imply that not even all the water 
behaves as a fluid. It is clear that any thickness of sediment could be obtained, at 
least within the usual range of values, by varying the assumed ratio of crustal 
thickness to fluid depth. 

Most authors have assumed, for purposes of computation, that the thicknesses 
of crust and fluid layer are equal. Curve fitting has then consisted of selecting the 
most favourable value of this thickness. The assumption is at variance both with 
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Fic. 14.—An example of the indeterminacy of Rayleigh wave solutions 
for oceanic paths. Dispersion from Prince Edward Island to Mirny 
(cf. Figure 8b), with three different interpretations. 
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the principle of isostasy (which seems to hold good in oceanic regions except for 
local anomalies), and with the evidence of seismic refraction results as summarized 
by Hill (1957). The case of a crust twice as thick as the fluid layer has occasionally 
been considered (see Press & Ewing 1955), but at the same time the very low 
value 5-5km/s has been adopted for the compressional velocity in the crust, in 
place of the usual 6-gkm/s. This has obscured the effect of varying Hi/Ho. The 
refraction results actually show that the average compressional velocity in the 
oceanic crust is very nearly the same for a thickness of 10km as for 5 km. 

We thus see that Rayleigh wave dispersion does not, by itself, provide a reliable 
measure of crustal thickness in oceanic regions, much less of the thickness of fluid 
sediments on the ocean floor. The high degree of uniformity that has been attri- 
buted to the oceans may be illusory in as much as it depends largely on Rayleigh 
wave evidence. Indeed the growing amount of data from seismic refraction work 
at sea, and also the evidence of Love waves, suggests that thicknesses of 1okm 
or more are not uncommon. We have found that this is the case in the oceanic 
regions adjacent to Antarctica. A contrary view of the major oceans regards the 
Pacific as having an essentially different crustal structure from the other oceans. 
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Since the Pacific is appreciably the deepest ocean we must expect it to have the 
least average thickness of crust, but the difference is not such as would suggest 
any uniqueness of structure. 


7. Structure of Antarctica 


It is concluded that eastern Antarctica as a whole is continental, for the present 
study has given a good coverage of this large area, and a crustal thickness of about 
35km has been consistently obtained. By contrast, Marie Byrd Land has an 
average crustal thickness of about 25km and cannot be regarded as truly conti- 
nental in structure. The data do not extend to the remainder of western Ant- 
arctica, nor has it been possible to study the important transition zone which runs 
from the Ross Sea to the Weddell Sea. 

These findings are evidently related to the known geological contrast between 
eastern and western Antarctica. The result for the western sector, however, re- 
quires further elaboration, for the average crustal thickness of 25km could be 
made up in a variety of ways. Beneath the high mountain ranges of Marie Byrd 
Land one would expect to find crustal roots giving locally a total thickness of at 
least 40km. If this is so there must be other zones where the crust is much thinner 
than 25 km, implying that the rock surface lies well below sea-level. The offshore 
depth contour at 1000 fathoms generally lies much further from the coastline in 
western than in eastern Antarctica; this would account for some small fraction of 
the difference in average crustal thickness. 

The picture of Marie Byrd Land as an archipelago made up of mountainous 
island chains finds support in recent seismic and gravity observations. A detailed 
survey of ice thickness has revealed extensive regions where the rock surface would 
undoubtedly be below sea-level if the ice were to melt. For example, along a 
traverse running close to latitude 77°S, between longitudes 87°W and 113° W, 
the rock surface was located consistently at between 1 km and 2 km below sea-level 
for a distance of some 400km (Bentley & Ostenso 1959). 

The conclusion that eastern Antarctica as a whole is continental and that the 
smaller average thickness of crust is confined to western Antarctica differs substan- 
tially from the conclusion of Press & Dewart (1959). For paths in all directions 
from Wilkes Station except to the southeast these authors obtained an average 
crustal thickness about as small as that indicated for Marie Byrd Land in the 
present study, and concluded that “at most only three-fourths of the Antarctic ice 
sheet is underlain by continent, the remaining area being oceanic in structure’’. 
The two exceptional paths, which passed close to Hallett Station, indicated a fully 
continental thickness. Six paths included approximately equal distances across 
eastern and western Antarctica, and traversed the Ross Sea-Weddell Sea zone. If 
it happens that this zone includes a wide belt of very thin crust, the average thick- 
ness obtained by Press and Dewart may prove to be consistent with the present 
results. ‘Three further paths crossed eastern but not western Antarctica, and for 
these it is not easy to find agreement unless there is a large region of oceanic crust 
south-west of Wilkes Station, which seems unlikely since Wilkes itself is situated 
on rock and the value of gravity there is close to normal. 
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Focal Mechanisms of some Earthquakes of the 
year 1950 


A. R. Ritsema 
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1. Data and Methods 


The data used in this study are compressions and dilatations of P and PKP 
waves to be found in the International Seismological Summary of the year 1950, 
supplemented with some initial motion directions of these waves asked for and 
contributed by seismic stations situated in the zones relevant to the determination 


(Table 1). 


Table 1 


Additional data not reported in ISS 
Earthquake 
No. Data 
(see Table 2) 
de Bilt C; Honolulu C; Hyderabad D. 
de Bilt D; La Paz (C’); New Delhi D; Scoresby Sund D. 
de Bilt D; Djakarta C; Honolulu C; La Paz C; New Delhi D. 
de Bilt D; La Paz D; La Plata D; San Juan C. 
La Paz C; 
de Bilt D; Djakarta D; Honolulu D; New Delhi D; Ottawa C. 
de Bilt C; 
College D; Resolute Bay C; Sverdlovsk D. 
Christchurch C; de Bilt D’; Djakarta D; Honolulu C; Wellington C. 
de Bilt D; Djakarta C’; Honolulu C. 
College C; de Bilt C; La Plata (C); San Juan D; Scoresby Sund C. 
de Bilt C’; Djakarta C; La Paz (C); Wellington C. 
Bogota D; de Bilt C; Guadalayara D; Honolulu (D); La Paz D; 
Scorseby Sund D. 
Bogota C; de Bilt D; Honolulu C; Scoresby Sund C. 
Christchurch D; de Bilt D’; Honolulu C; Hyderabad D; Riverview D; 
Sendai D; Wellington D. 


I 
2 
3 
4 
5 
6 
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Remarks: C : compressional P wave 
: dilatational P wave 
”: compressional PKP wave 
‘: dilatational PKP wave 
): inconsistent with the solution reached. 


The amplitude of the longitudinal wave radiated as a result of a sudden motion 
along a flat fault plane in the focus of an earthquake is proportional to 


sin 26 . sin ¢, 
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where @ and ¢ are spherical coordinates measured relative to the direction of the 
fault motion and the position of the fault plane respectively (Hodgson 1957). 
There are thus two quadrants with compressional onsets of the longitudinal wave 
interchanging with two quadrants with dilatational onsets. In the fault plane 
(¢@ = ©) and in the plane perpendicular to the direction of fault motion (@ = go°) 
the amplitude of P is reduced to zero; these are the so-called nodal planes for longi- 
tudinal waves. 

Assuming now that the initial character of the longitudinal wave remains un- 
changed on the pathway from the focus to the recording station, it is possible to 
determine the position of the two nodal planes in the focus with the help of the 
compression or dilatation data of many seismic stations all over the world. The 
problem is reduced to the determination of the angle 7, with the downward vertical 
in which the wave left the focus to emerge in the recording station at the known 
epicentral distance A. This (i,, A)-function has been determined with the help of 
the travel-time tables of Jeffreys and Bullen ((t, A)-function) and the velocity 
distribution of longitudinal waves in the upper mantle derived by Jeffreys ((vr, r)- 
function), by means of 


where R is the Earth’s radius, r is the radial distance from the focus to the centre 
of the Earth, v; is the velocity of longitudinal waves at the depth of focus, and 
dA/dt is the velocity of the P wave front along the Earth’s surface (Ritsema, 
1958). 

The compression and dilatation data are now plotted in a stereographic 
diagram. The tangential co-ordinate in the diagram is the azimuth of the re- 
cording station as seen from the epicentre, and the radial coordinate is taken in 











proportion to the value of tan $7. The compressions and dilatations are separated 
by two nodal lines that both follow the course of a meridional line in Wulff ’s 
stereographic net and that in addition stand perpendicular to each other. The 
position of the nodal planes in the focus represented by these lines can directly be 
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read in the stereographic net. One of these nodal planes is the fault plane and the 
motion direction in the plane is given by the pole of the auxiliary plane. 


2. Presentation of the results 


The solutions, presented in Table 2, are given in the form of the positions, 
(azimuth and plunge) of the A-, B- and C-axes and of the earthquake type. 

The fault motion direction is called the A-axis, the line in the fault plane per- 
pendicular to the A-axis is called the B-axis, and the normal to the fault plane the 
C-axis (Figure 1). The B-axis is determined solely by the longitudinal wave data, 
it is the line of intersection of the two nodal planes. The positions of the A- and 
C-axes—the poles of the two nodal planes—are also determined in this way. 
These two axes, however, cannot be distinguished mutually without an investiga- 
tion of the initial motions of the transverse waves in many directions around the 
focus. Since S waves have not been considered, the solutions in Table 2 are am- 
biguous with respect to the A- and C-axes. For each earthquake both possible 
solutions are given (a and b respectively). 

The earthquake type is given in the form of a combination of the symbols 
P and T with R and L according to the accompanying scheme (Figure 2). The 
first of the two symbols noted gives the dominating characteristic of the earth- 
quake. 

Earthquakes of mainly P and mainly T type are characterized by a B-axis 
plunging at a smaller angle than 45°. The B-axis is the direction in which the 
second greatest stress component (gc) is acting. ‘That means that the so-called 
plane of action, formed by the A- and C-axes, in which plane also the two principal 
stresses (oc; and og) are situated, dips steeper than 45°. Assuming that the angle 
with the principal stresses in which the fault motion takes place is approximately 
45°, type P is defined in such a way that the maximal stress component (0) is 
more nearly horizontal than the minimal stress component (a3). In type T the 
minimal stress component (3) plunges at a smaller angle than the maximal stress 
component (0). Motions that are characterized by P and T type stress configura- 
tions are along reverse and normal faults respectively. 

In the case that the B-axis plunges steeper than 45°, the principal stress compo- 
nents (o; and o3) are situated in a plane dipping at a smaller angle than 45°. ‘Trans- 
current fault motions are characterized by such a stress configuration. The earth- 
quake type is given as R or L respectively according to the dextral or sinistral 
character of the resulting fault motion. 

Pure normal, pure reverse and pure transcurrent fault motions only rarely 
occur. Mostly there is a smaller or greater transcurrent component in mainly 
normal and mainly reverse fault motions; and also a smaller or greater normal or 
reverse component in the mainly transcurrent fault motions. The characterization 
of this minor component is given as a second symbol in the earthquake type. 


3. Conclusions 


The predominance of transcurrent fault motions in earthquake foci over reverse 
and normal fault motions is a confirmation of older results (Hodgson 1957). 
That the A- and C-axes show a preference for directions that are about perpendi- 
cular to and parallel with the seismic zones in which the earthquakes take place 
also is in accordance with known facts. 
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The principal stresses are denoted by o1, o2, and os, where o1 > G2 > Gs. 
The A-axis is the fault motion direction, the B-axis is the line in the fault 


plane perpendicular to the A-axis, and the C-axis is the normal to the 
fault plane. 


The B-axis is parallel to oe. 
The plane of A- and C-axes is identical with the plane of o1 and os. 


The fault types are given as an illustration of the different stress configura- 
tions rather than as an equivalent definition for the earthquake types. 
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A Theory of Polar Geomagnetic Storms 


J. H. Piddington 


(Received 1960 April 20) 


Summary 


This is the third of a series of papers outlining a hydromagnetic 
theory of geomagnetic storms. The first two dealt with the storm-time 
variations Dst in terms of isotropically propagating hydromagnetic 
compression (first phase) and rarefaction (main phase) of the geo- 
magnetic field and its plasma. 

The three phases of the longitude-dependent disturbance Ds, the 
first phase or Sc, the main phase and the preliminary reverse are all 
accounted for in terms of hydromagnetic waves whose Poynting fluxes 
are directed along the magnetic lines of force. These disturbances 
originate from “electromagnetic friction” between the solar wind and 
geomagnetic lines of force near 06 and 18h local time. These lines 
of force are bent out of their meridian planes and cause two pairs of 
“twist” waves to reach the Earth, being focused into relatively confined 
regions at high latitudes. These twists are stable or force-free and 
result in space-charge accumulations in the lower ionosphere where 
the resultant Hall current accounts for the observed disturbances. 

During the main phase the twists are maintained by a deformation 
of the geomagnetic field which was described in the first two papers 
(the geomagnetic “‘tail”). Thus the tail theory accounts for both Dst 
and Ds main phase. It is difficult to see how any ring-current theory 
or other theory dependent on trapped corpuscular radiation could 
account for the asymmetric features of a storm. 

The mechanism also provides the beginning of a theory of aurorae. 


1. Introduction 





It is generally agreed that geomagnetic storms and associated effects result 
from the passage past the Earth of neutral, ionized gas ejected from the Sun. 
Chapman & Ferraro (1931) showed that a cloud of such gas approaching the 
Earth would have electric currents induced in its surface layers and they explained 
the world-wide sudden increase in the horizontal component of the geomagnetic 
field as a direct result of these currents. 

The discovery that the outer geomagnetic field encloses a highly conducting 
plasma has necessitated some modification of this theory. The currents at the 
surface of the cloud cannot affect the field directly (that is according to the Biot- 
Savart law) at the Earth. They can, and do, initiate a hydromagnetic disturbance 
comprising a general movement and distortion of the field and of the plasma into 
which the field is “frozen”. This disturbance travels into the vicinity of the 
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Earth with a speed of ~107—10®cms~! until it envelopes the lower ionosphere. 
The geomagnetic disturbance then diffuses through the ionospheric E-region 
as through a solid conductor according to ordinary “‘skin-effect”’ theory (Pidding- 
ton 1959, called P.1). 

Apart from this transmission theory the important difference between the 
Chapman-Ferraro theory and the hydromagnetic theory is that solar gas cannot 
flow freely into the geomagnetic field along lines of force which cut the Earth at 
high latitudes (P.1). These “thorns” of gas are rigorously excluded by violent 
self-inductance effects and all the lines of force are bent (P.1, Figure 2), rather 
than cut (Chapman & Ferraro, Figure 4). 

A second difference between the two theories follows from this: if the solar gas 
may only diffuse into the geomagnetic field then the initiation of the mass motion* 
required by a large-scale ring current seems improbable. Chapman & Ferraro’s 
(1933) theory of the main phase of the Dst disturbance and Martyn’s (1951) auroral 
zone theory encounter this difficulty ; also it is likely that the fast westward drift at 
~108cms~! required by this theory would have been noted during satellite 
experiments. An alternative mechanism involving a geomagnetic “‘tail” on the 
dark side of the Earth has been proposed (Piddington 1960, called P.2). 

These theories have all dealt with the world-wide storm-time variations Dst, 
averaged in longitude and comprising an increase in the horizontal component of 
the geomagnetic field during the initial phase (SC) and a decrease during the main 
phase. When the Dst variations have been eliminated the remaining disturbance 
(Ds) is found to consist of three separate phases, each highly dependent on longi- 
tude as well as latitude and each concentrated towards the polar regions. It is with 
these complex systems of disturbances, which may be referred to broadly as 
“polar magnetic storms’’, that the present paper is concerned. 


2. Polar magnetic storms 


The total geomagnetic disturbance during a sudden commencement storm is 
represented by a function D of the longitude Ag referred to the direction of the Sun 
and the time tst reckoned from the beginning of the storm. Thus, as Obayashi & 
Jacobs (1957) have written: 


D(ds; tst) = Dst(tst)+ Ds(Ag; tst). 


The longitude dependent part of the storm Ds takes three distinct forms during 
three phases of the storm as follows: 


(a) During the initial phase or sudden-commencement (SC) phase the field 
has a form referred to here as Ds(SC). 

(6) During the main phase disturbance field changes to Ds(MP). 

(c) Sometimes, prior to SC, there is a preliminary reverse disturbance (SC*) 
starting a minute or so before SC. This is referred to here as Ds(PR). 


The longitude dependent part of magnetic storms which have no sudden 
commencement generally have a form similar to Ds(MP) and the theory of their 
origin is similar to that given below. 


* This motion is not merely the differential drift which constitutes the ring current. In addi- 
tion the whole mass of plasma within the current zone must spin around the Earth with a speed of 
order 10°cms~!. This mass motion and the centrifugal force resulting from it is the basic require- 
ment of the main phase of a storm; without it the current itself would rapidly die away. 
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The averaged atmospheric current systems which could account for these 
disturbances are complex; examples calculated and plotted by Obayashi & Jacobs 
(1957) are reproduced in Figure 1. These diagrams indicate that in polar regions 
the horizontal intensity could either increase or decrease during any of the periods, 
depending on Ag and on the latitude. The whole of each of the three equivalent 
current systems is needed for a satisfactory description of the storm phenomenon, 


although it may be that Ds(PR) and even Ds(SC) are absent from particular 
storms. 


Fic. 1.—Atmospheric current systems which could account for the 

three phases of the Ds geomagnetic disturbance field (after Obayashi & 

Jacobs). (a) The SC phase, (b) the main phase, and (c) the preliminary 
reverse phase. 


An examination of Figure 1 suggests that the main features of the Ds variations 
are as follows: 


(a) The disturbances comprise sets of clockwise and anticlockwise current loops 
situated mainly in and above the auroral zones with one system of loops in the 
morning and one in the evening hemispheres. 

(6) The disturbance vector changes rapidly both with longitude and latitude, 
within distances less than 1000km. This is remarkably small scale considering 
the likely dimensions of the original disturbance (~ 30 000 km). 

(c) The Ds(MP) currents are generally similar to the Ds(.SC) currents except 
for greater concentration in the auroral zone and also for two weak lower-latitude 
current loops. 

(d) The Ds(PR) currents are almost perfectly reversed and weakened Ds(SC) 
currents (see Figures 1a and 1c). 

(e) ‘The Ds(PR) appears to be the initial part of the same disturbance which 
causes Ds(SC); when Ds(PR) is present then Ds(SC) is delayed for the duration 
of the reverse impulse (Gerard 1959). 

(f) The centre of the morning disturbance occurs at about 05 hours local time 
or earlier; that of the evening disturbance at 17 hours or earlier. 

(g) The intensity of the disturbance is about four times that of the Dst field. 


3. The location and immediate cause of the Ds current systems 


The equivalent (atmospheric) current systems of Figure 1 are not necessarily 
the systems actually responsible for the Ds field—they are merely a handy way of 
displaying this disturbance. The location of the actual Ds currents is the first 
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step in formulating a theory of their origin, but first it is instructive to review the 
corresponding problem for the Dst field. 

There can be little doubt that the first phase of Dst results from a compression 
of the plasma and its frozen-in geomagnetic field. The compression wave travels 
towards the Earth causing a general increase in the horizontal component of the 
geomagnetic field. On reaching the ionosphere the wave sets up currents which 
account for the initiation of the observed geomagnetic disturbance. However, 
because of the low integrated conductivity of the ionosphere the electromagnetic 
wave front soon penetrates this region and establishes itself as a current sheet in 
the body of the Earth. Thus the Dst field, while originating in ionospheric current, 
is maintained by two other current systems*, one in the Earth and one adjacent to 
the solar wind (P.2). 

Similar considerations apply to the main phase of Dst when the horizontal 
component of the field decreases everywhere due to an inward propagating rarefac- 
tion wave. The important point is that in either case ionospheric currents are 
driven by an induction electric field which only lasts while the requisite magnetic 
flux is entering or leaving the ionospheric “‘shell’”’ and then ceases. 

The Ds current system is much more complex but, as before, no currents above 
about 1000km may affect the field at the Earth’s surface directly (according to the 
Biot-Savart law). This result is consistent with the tentative conclusions of 
Chapman & Bartels (1940, Chapter 22) based on harmonic analysis of the surface 
field. More distant current systems may, however, initiate hydromagnetic motions 
which result in currents in the ionosphere or Earth or both. The alternative, that 
ionospheric currents are due to winds there, appears most unlikely from the 
discussion in the appendix. 

Let us then consider a possible origin of Ds in terms of a hydromagnetic dis- 
turbance incident on the lower ionosphere from above. The equivalent current 
systems indicate that its form must be quite different from that of the Dst com- 
pression and rarefaction waves. To determine its form consider an elementary 
ionospheric current loop as shown in Figure 2, corresponding to one current loop 
around region A of Figure rb. An explanation of this elementary loop could be 
extended to cover all the main features of all three phases of the Ds field. 

The clue to the origin of the current loop is the nature of the electric field 
required to drive the current. If the medium concerned were isotropically conduct- 
ing, this would be necessarily an induction field E; in the direction of the current 


and given by 


curl Ej = — div E; = 0, (1) 


dt’ 
where H is the magnetic vector. The hydromagnetic disturbance which could 
cause such a field comprises the expansion of a cylindrically shaped bundle of 
magnetic lines of force. It is shown in the appendix that such disturbances could 
not account for the main Ds features and that the actual current system must flow 
in the ionosphere rather than the Earth. 

In the ionosphere the conductivity is anisotropic, having components oo along 
the magnetic field, and o; (Pedersen) and o2 (Hall) perpendicular to the field. In 


* The statement in P.1 (Summary) that “the currents responsible for all geomagnetic dis- 
turbances must flow at levels below about 1000km” should be corrected to “the currents directly 
responsible...’’. ‘The Earth is certainly shielded from the direct effects of the current sheet at the 
surface of the solar wind. However, as already stated, the latter initiates the train of events leading 
to the observed disturbances and it is also needed to maintain any quasi-stable Dst field. 
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polar regions (and in Figure 2), the currents flow nearly perpendicular to the field, 
so that we are concerned with o; and o2 and, when Hall current is prevented by 
the development of space-charge, an “‘effective” conductivity ¢3 = 0; +(022/01) 
(Baker & Martyn 1953). An electric field without divergence (equation (i)) 
would cause the observed current loop* but as seen in the appendix such a field 
may not be invoked because it may only be applied for brief periods (correspond- 
ing to the passage of a hydromagnetic disturbance). 

A second type of electric field may also explain the current loop in Figure 2, 
provided this flows in an anisotropically conducting medium. This is a potential 
field Ep radially outwards and due to an accumulation of positive space-charge 
near O, having density g(e.m.u.). Then 


curlE, = 0, div Ey = 47qc?, (2) 


and since og > oj, the current is mainly Hall current. Since it flows in closed 
loops the maximum effective conductivity og is not developed. As well as Hall 


j 


Fic. 2.—Two possible ionospheric electric field systems which could 
explain the Ds equivalent current loops. 


current, Pedersen current o; Ey must flow radially. This in turn requires a current 
to flow downwards into the ionosphere in the interior of the loop (and in region A 
of Figure 1b) and upwards out of the ionosphere near the boundaries of the loops. 

The nature of the hydromagnetic disturbance needed to drive the Ds current 
system is now evident. It comprises a twist in the bundle of magnetic lines of force 
passing through the region between A and the outermost current loop of this 
particular (clockwise) system. This twist is identified with the downward flowing 
current j by the equation curl H = 47j. This disturbance differs in two important 
ways from the Dst disturbances. First it is substantially force-free because the 
current follows the lines of force; this means that the phenomenon is magneto- 
hydrostatic and stable. Second, the downward current cannot flow below the 
ionosphere so that the disturbance remains anchored there instead of descending 
into the Earth. 

In the following section it is shown that there are two distinct types of hydro- 
magnetic wave which may propagate above the ionosphere. One of these is 


* The appropriate conductivity would be ¢3 because the Hall current o2H a (E,/H) tends to 
flow radially and so is prevented by the development of a space-charge field. 
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responsible for the Dst disturbances and the second for the Ds disturbances. The 
latter has the distinction that it may only propagate along magnetic lines of force 
so that it tends to concentrate its effects near the poles, as observed for the Ds field. 


4. The two types of wave responsible for the transport of the Dst 
and the Ds disturbances 


There are two types of weak hydromagnetic wave*, having quite different 
properties. The differences are revealed by the equation for the mechanical force 
exerted by the distorted magnetic field on the plasma: 


F=jaH= “(curl H) AH 
47 (3) 


(He I 
= —grad (—) — : V)H. 


Of the two stress terms the first represents a gradient in the magnetic pressure 
and permits propagation as for a sound wave; the second denotes ‘‘magnetic stiff- 
ness” and permits propagation as for a transverse wave along a stretched string. 

One type of hydromagnetic wave has the perturbation magnetic field h and 
and the gas velocity v in the plane containing H and the wave normal n. This wave 
uses either or both of the above stresses and always propagates with velocity 
V = H(4zp)-+*, where p is the plasma mass density. The energy flow (direction 
of the Poynting vector) is along n. 

The second type of wave has h and v perpendicular to the H-m plane and 
this removes the stress term grad(H?/87) so that the wave is purely transverse 
and propagates with velocity V cosy where # is the angle between H and n. 
The Poynting flux is along H so that the energy flow is directed along the main 
field rather than along the wave normal. 

The velocities of these two wave types, called Et or “isotropic” and O or 
“twist” waves respectively are shown in Figure 3 (reproduced from Piddington 
1955). The most important difference between them, as far as geomagnetic dis- 
turbances are concerned, is the direction of energy flow. 

In the geomagnetic meridian plane containing the Sun, the solar wind must 
distort the field always so that h also lies in this plane. A disturbance moving 
towards the Earth has n in the same plane and is therefore an isotropic wave. The 
result is a more or less uniform change in magnetic energy around the Earth, an 
increase during the initial phase and a decrease during the main phase. This is 
the proposed extension and modification (P.1, P.2z) of the Chapman & Ferraro 
theory to account for the Dst field. 

In meridian planes within a few hours of sunrise and sunset the solar wind 
moves mainly at right angles to the planes containing the magnetic field lines. 
There is a tendency here, as noted by Chapman & Ferraro and Martyn (1951), 
for the solar wind to blow through the outermost regions of the field. This is 
permitted by the setting up of an electric space-charge field which opposed the 
induction field v A H caused by the passage of the plasma through the geomagnetic 


* A third wave which may be regarded as a modified sound wave is probably not important 
in the exosphere because of the low gas pressure. 

t The symbols E and O stand for extraordinary and ordinary; at high frequencies these waves 
become the two well-known radio waves in a magneto-ionic medium. 
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Fic. 3.—Wave fronts of the E or isotropic and O or twist hydromagnetic 
waves originating in the region near O. 


ahs 





of 





solar_wind v 


Fic. 4.—A schematic diagram of the Earth’s magnetic dipole field 

viewed from a great distance above the north geomagnetic pole. Near 

the geomagnetic equatorial plane the solar wind of velocity v flows in 

the regions beyond the dashed lines (one marked RR). Electric space- 

charge fields Ep are set up and the magnetic lines of force bent as shown 
by the dot-dash lines. 





A theory of polar geomagnetic storms 321 


field. In the morning hemisphere solar protons are forced inwards towards the 
Earth and electrons outwards, thus creating an outwardly directed space-charge 
field; an inward field develops in the evening hemisphere. 

The situation is shown schematically in Figure 4 where the Earth is viewed 
from a long way above the north geomagnetic pole. The magnetic lines of force 
are all in meridian planes passing through the point of observation and so appear 
as straight lines, longer ones cutting the Earth at higher latitudes as shown. The 
inner boundaries of the solar wind, in regions not far from the equatorial plane, 
are represented by dashed lines on which electric space-charge is shown. Space- 
charge of opposite sign accumulates beyond the boundaries of the solar wind so 
that the potential electric fields shown as E, are set up according to equation (2). 
The Lorentz forces responsible for setting up Ep are very powerful but as E, tends 
to neutralize E; = va H, they are balanced and the situation tends to be stable 
without further flow of current or dissipation of energy of the solar wind. Situa- 
tions where the dissipation is zero have been examined by Chapman & Ferraro 
(1931) for an infinite slab and by Ferraro (1937) for the non-uniform rotation of 
a magnetic star (the case of isorotation). 

In the present situation the electrically polarized slabs of gas are not infinite 
in extent and substantial space-charge must leak away all around their boundaries. 
The current density in the slab (j) required to replace this loss flows in the direc- 
tion of E; and exerts a force on the solar wind given by equation (3) and opposite 
to v. Thus the current leakage results in a frictional drag on the solar wind. Addi- 
tional drag may result from the formation of small space-charge clouds in the 
region of interaction (P.2) and by hydromagnetic instability (Parker 1959). 

For the purpose of the present discussion the important result of the drag on 


the solar wind is that it causes a large-scale distortion of the geomagnetic field 
which propagates along the magnetic lines of force until it reaches the Earth as a 
twist hydromagnetic wave. This is discussed in the following section. 


5. The electromagnetic field components of the Ds disturbance 


The electric current systems in and around the electrically polarized slabs of 
solar wind are complex, particularly near the outer limits of the field (as at P, 
Figure 4). Nearer the centres of the slabs however (Q) the current should be pre- 
dominantly in the direction vAH of the Lorentz force, even though this is 
opposed by a weaker space-charge electric field*. ‘This current system, which 
for weak perturbation is everywhere nearly perpendicular to Hp and to the vector 
Vv, initiates the sequence of events leading to Ds. 

In Figure 5 the geomagnetic field is shown schematically, as seen from a direc- 
tion a few degrees north of the Sunt. Two undisturbed lines of force are shown 
(full lines) and this field is referred to as Ho. The solar wind passes through the 
outer parts of this field and the wind current system jy» is represented by the full 
lines at the left of the figure. The perturbation magnetic field due to jy is hy, having 
the form of loops with hy» clockwise when looking towards the Earth. One such 
loop is shown, seen by an observer who has moved from his original position 


* The question of Hall current does not arise because the electrons and positive ions do not 
spiral. If the opposed fields were equal, then the ions would follow straight paths. To replace 
boundary leakages the positive ions drift to the right across the lines of force in the direction 
v AH. A possible secondary effect of Hall current is discussed in Section 9 below. 

t So that the distorted lines of force may be shown in perspective. 
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towards the left-hand side of the diagram. This current system has no direct 
effect near the Earth because of shielding. 

Near the geomagnetic equatorial plane the field h,, opposes Hp in front of the 
plane of the paper and reinforces it behind that plane and so constitutes a move- 
ment of lines of force away from the Sun. In the (meridian) plane of Figure 5 
and north of the equatorial plane h,, is towards the Sun, south it is away from 
the Sun. Clearly the combined field H = Ho+h, is as shown by the dot-dash 
lines in Figure 4, the lines of force being on surfaces formed by bending the outer 
parts of meridian planes away from the Sun. Such a bend in the lines of force 
must propagate as a hydromagnetic disturbance by virtue of the force term 
(H* /)H of equation (3). At the bend the perturbation field is perpendicular to 





Fic. 5.—A hydromagnetic twist caused by the solar wind dragging on the 
geomagnetic field. To provide perspective the Earth is viewed from a 
point north of the Sun. 


the meridian plane which contains the main field Hp. Any wave which is to reach 
the vicinity of the Earth and cause observable geomagnetic effects must have its 
normal in that plane and so must propagate as an O wave with its Poynting flux 
along Hy (Section 4 above). This explains why the observed Ds disturbance is 
mainly polar. 

Now let us consider the hydromagnetic O wave when its front has advanced to 
the position WF shown by dot-dashed lines in Figure 5. It must have separated 
into two parts, one following the lines of force to the northern and one to the 
southern hemisphere. Just as for the Dst wave front, the passage of the Ds wave 
front must cause changes in the geomagnetic field and its plasma as it passes 
through a particular region. These changes will remain after its passage so that 
before examining conditions inside WF it is desirable to consider those between 
WF and the solar wind. 

Consider the line RR of Figure 4 representing the inner boundary of the solar 
wind. Outside this boundary the current j» is mainly across the magnetic lines of 
force. Inside the boundary, however, the plasma is at rest and the current j consti- 
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tutes a leakage of space charge from the wind region. Because of the much higher 
conductivity along the magnetic field the current flows mainly in that direction 
as shown in Figure 5. This current system not only accounts for leakage from RR 
and hence the wind current system itself but also, as shown below, is a necessary 
adjunct to the wave front itself. The magnetic field h caused by j is similar to 
h,, but it combines with Hp in a very different manner causing the lines of force 
to assume spiral forms. In other words j introduces a twist into the bundle of 
lines of force in which it flows. In the northern hemisphere and looking along Ho 
it is a right-handed twist. 

The electromagnetic field within the wave front itself is examined in the 
following sections. It is evident, however, that as WF defines the lower boundary 
of the twist, current must flow down the magnetic lines of force into its central 
region and then out across the lines of force and in the plane WF. Some must, 
escape into interplanetary space as shown by the vector J; some must flow in other 
directions and so into the geomagnetic field. Variations in plasma conductivity 
and extent of the current paths make it likely that the J current will vary around 
the wave front so that the twist is asymmetric. 

The three current systems discussed above and their magnetic fields and 
mechanical forces may be described in a more quantitative manner by use of the 
equations: 


curl hy = 47jw, curl h = 47j. (4) 


Consider a right-hand system of axes with Ox tangential to Hp at the origin and 
Oy in the direction v; the directions of Ox and Oz change as O moves about in 
the geomagnetic field. In the region penetrated by the solar wind we have 
jw ~ ©, O,—Jwz, so that 

Chwz Chwy 


by oie = 47) wz, (5) 


which indicates that the lines of force are bent like the dashed line in Figure 5. 
There is a retarding force of density fy = —jwzHo on the wind and an opposite 
stress on the magnetic field. Near the equatorial plane the important perturbation 
is hw, so that 


hw 
fy ~ = = . (6) 

47 oy 
Integrating in the Oy direction, the total force per cm? opposing the solar stream 
is Hohwz/47. Far from the Earth the dipole field is too weak to oppose the solar 
stream and hy»z must be small. It should increase to a maximum in the region 
where the perturbation field is approximately equal to the dipole field, that is 
where the lines of force are substantially bent. Here the Awz ~ Ho ~ hy» and the 
total force hy,2/47 must be equated to the rate of loss of momentum of the solar 

stream: 

hy? = 8e8BNM2?. (7) 
Here N, M and ware the ion density, mass and velocity and 8 (< 1) the proportion 


of the total momentum lost. It is shown below that the relative intensities of Ds 
and Dst suggest a value of B < 0-4. 


Between the solar wind and WF the important current is along the field, so 
that 


(8) 
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which describes a twisted bundle of lines. The force on the plasma is jh H = 0 
so that the field is force-free. Since there are no electromagnetic forces on the 
plasma (which would need to be balanced by inertial forces or pressure gradients) 


the situation should be stable. Conditions inside WF are discussed in the following 
section. 


6. The propagation of a hydromagnetic twist 


The Ds twists are probably asymmetrical but their important features should 
be revealed by analysing a circularly symmetrical twist shown schematically in 
Figure 6, which is now described. 














Fic. 6.—A simplified (cylindrical) hydromagnetic twist propagating down 
the magnetic field Hz. 


In a right-handed system of cylindrical co-ordinates (r, 9, z) assume a uniform 
magnetic field and ion plasma contained within a cylinder r = R. The field is 
parallel to Oz and may be assumed substantially frozen into the plasma which, in 
turn, is contained within an insulating cylinder. A twist is now introduced into 
part of the bundle of lines of force and the propagation of the disturbance is 
examined. 

Assume that in the region z < o there is a perturbation magnetic field 
h = 0,h,{r),o which is associated with a current of density j = 0, 0, j-(r). 
Maxwell’s equation relating these parameters is then 


10 ; 
—-—(rhg) = 4nj2. (9) 
r or 


There must be a second current system /J;,(r), being the total radial current density 
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in the surface of discontinuity z = o. The equation divj = o must be satisfied 


to a very high degree of approximation and integrating through the region of 
discontinuity we have 


1@ ; 
-—(rJr)—jz = 0. (10) 
r or 


Combining these equations we find 
hg = 47Jr. (11) 


The propagation of the boundary of this twist may be studied by transforming 
to axes in which the magnetic field does not change with time, so that the wave 
front is at rest and the plasma enters it with velocity —v, and leaves with velocity 
(0, —vs, —vz), as shown in Figure 6. The force-momentum equation in the 6 
direction is then 

pvz9 = JrH:, (12) 


where p is the mass density of the plasma. Also, if the gas is constrained to move 
along the magnetic lines of force, then 


h 
U = vs. (13) 
H, 


Combining equations (11), (12) and (13), 


vz = HA4n7p)+, (14) 


so that the front of the twist advances through the plasma with the usual Alfvén 
velocity. 

This result shows that twist waves have velocities which are independent of 
their wavelength and strength (provided they remain weak). The medium is non- 
dispersive and the disturbance advances as a single front with conditions changing 
more or less abruptly from one side to the other. The propagation of an asym- 
metric wave, that is a wave in which A, varies with 6, should be similar to that 
analysed above although its effects on arriving in the ionosphere will be different. 


7- Electric fields and “slipping” 


In the previous section propagation of the twist was discussed in terms of the 
momentum change of the plasma and the consequent electric current and magnetic 
field gradient. It is also necessary to consider the electric field responsible for 
Jr. If the plasma were isotropically conducting this would simply be the induction 
field caused by the flow of plasma across the perturbation field h, in the transition 
region (Figure 6): 


E; = vAH = (v,hy, 0, 0). (15) 


This field drives current out of the twist cylinder where, to ensure its further 
progress (to complete the current loop), positive space charge will accumulate, 
giving rise to a potential field Ep. 

This secondary field has several interesting effects: in the first place it per- 
forms the function for which it was created, namely of driving current up the 
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lines of force surrounding the twist thus limiting the extent of the twist and allow- 
ing its current to escape. Within the wave front Ep, is radially inwards, thus op- 
posing E;. The difference between the two is related to the total current J; inte- 
grated through the front, by Ohm’s law: 


Jr = = (vetg—Ep), where = = Jods, (16) 


o being the ordinary conductivity and = the conductivity integrated through the 
slip region. 

If the external circuit has high resistance then Ep» will largely cancel E; as in a 
homopolar generator with little or no load. The transition region then tends to be 
extensive in order that 2 may be large enough to provide the required current J; 
with a much reduced net electric field. Since it is the front of the disturbance 
region which travels at the velocity given by equation (14), the main part of the 
twist tends to be delayed behind this front. This delay and the prior arrival of a 
space-charge electric field is shown to be the cause of the preliminary reverse 
impulse Ds(PR). 

Now consider a plasma which is not only anisotropically conducting but in 
which the ion collision distance and gyro-radius may exceed the dimensions of the 
disturbance. The concept of conductivity may no longer be applied; rather must 
the paths of individual ions be studied. 

In axes in which the disturbance is at rest the plasma flows with velocity — wv, 
into the twist. The electrons may be regarded as firmly attached to the magnetic 
lines of force. The heavy ions, having a velocity component vzh,/H, perpendicular 
to the lines of force are deflected outwards and if space-charge could be neglected 
they would be deflected inwards again and would spiral down the lines of force. 
Thus the whole mass of the gas would be deflected along the lines of force as 
required. 

However, as the ions start to move out from the plasma* they create space 
charge. If they moved out one gyro-radius v,h,M/H *e, where e and M are their 
charge and mass, then space charge per cm? area of cylinder wall would be 


NMv-h, UN, 


a ete. P 
a a (7) 





The corresponding radial electric field is found by integrating equation (2): 
Ep = (Epr(r), 0,0) = —47g0? 
—c (18) 


ve . Ugh. 
The induction field vzh, is responsible for the creation of the reverse (space- 
charge) field. The latter may approach the former in magnitude but cannot exceed 
it; yet if the ions were allowed to spiral freely then Ep, as given by equation (18), 
would exceed vzhy by the very large factor c2/V2. ‘This means that the plasma 
must continue to move nearly parallel to the axis of the cylinder; in other words 
the space-charge provides slip between the field and plasma as in the case of the 
isotropically conducting fluid. Since the space-charge field must cause external 
currents to flow there will also be some twist as well as the slip. 


* The enclosing insulating cylinder might here be assumed to be replaced by plasma at a 
pressure sufficient to contain the magnetic field. 
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Now consider the situation when the disturbance collides with a rigid semi- 
conducting sheet behind which is an insulating medium. The sheet is meant to 
represent the ionosphere which, as far as the propagation of geomagnetic dis- 
turbances is concerned, may be considered a rigid sheet (P.1). 


8. The principal ionospheric current system 

Although the preliminary reverse impulse, when it occurs, precedes the 
Ds(SC) and Ds(MP) disturbances, its discussion will be delayed to the following 
section. 

Consider the situation when the twist has been established right down to the 
lower ionosphere. There is no longer any slip, the current jz flows into the iono- 
sphere and the geomagnetic field lines are twisted (hg) down to the same level. 
The situation is magneto-hydrostatic, being stable because the current flows very 
nearly along the field so that the force on the plasma, ja H, is nearly zero, even 
though j is finite. This is the important difference between the Ds current system 
and the Dst system; the latter flows at a large angle to H and so moves the plasma 
in such a way as to transfer the current down through the ionosphere and into the 
Earth. 

The downward current system jz cannot continue past the ionosphere but must 
flow radially outwards from the cylinder and back up the magnetic lines of force 
which surround the cylinder. Thus a quasi-stable current is established in and 
above the ionosphere. However, as this current flows in a closed system of sheets* 
its magnetic field at the surface of the Earth is likely to be small. A more important 
feature of the radial ionospheric current is that a space-charge electric field must 
be established to drive it. The current flowing down from above causes positive 
space-charge to accumulate in the cylinder, giving a radial outward electric field. 
This field will drive a Hall current in a clockwise direction (viewed from above). 
Since the height-integrated ionospheric Hall conductivity is about twice the 
Pedersen conductivity (Baker & Martyn 1953, p. 291) the Hall current J, is about 
double J,. The magnetic effects of J, at the Earth are likely to be more than 
double those of /J;. 

It is this clockwise current system which is proposed as the cause of the first 
phase of a sudden-commencement geomagnetic storm in the morning hemisphere. 
A corresponding anticlockwise system, caused by an oppositely directed twist, 
would then account for the afternoon disturbance. 

In Figure 5 the outer line of force is shown as deflected out of its meridian 
plane away from the observer. If it is gripped firmly enough by the solar wind it 
will become frozen into the wind and be carried away to contribute to a geo- 
magnetic “tail”. This tail may explain the Dst main phase (P.1 and P.2) since it 
causes a general decrease in magnetic energy and magnetic pressure throughout 
the remaining geomagnetic field. It may also explain that part of the Ds field 
which continues after the solar gas cloud has passed the Earth. The magnetic 
lines of force will remain bent in much the same way as shown in Figure 5 except 
that the outermost part will be carried away entirely. 

In P.2 certain advantages of the magnetic tail theory of Dst over ring-current 
theories were discussed. However, the total main-phase disturbance comprises 
Dst+Ds, their separation being quite arbitrary; any satisfactory theory must 


* Comprising pairs of co-axial cylinders with current flowing down the inner cylinder and up 
the outer cylinder and connected via the ionosphere. 
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explain the total disturbance. Any ring current theory or other theory involving 
ions trapped within the geomagnetic field could only cause a circularly symmetrical 
disturbance such as Dst. Some separate mechanism, in which the direction of the 
wind vector v retained its significance after the wind ceased to blow, is needed to 
explain Ds. The magnetic tail theory explains the whole phenomenon. 

The equations of Section 5 may now be used to find the electric current for a 
given amount of twist of the magnetic field. For simplicity assume the pitch of 
the magnetic spiral to be independent of r so that 


hr) = rho, (19) 


where ho is independent of r. Then from equation (9) 


. to 
jis? (20) 
27 
so that jz is also independent of r. The ionospheric current is then found from 
equation (10) 


(21) 


During a moderate magnetic storm a total Hall current of ~2x 105A flows 
in a loop of radius ~ 1000km so that the current per cm strip J, ~ 2x 107-3 A 
or 2 x 10-4e,.m.u. If the integrated Hall conductivity is 1-4 x 10-8 (Baker & Martyn 
1953) then the electric field required is 1-4 x 104e.m.u. or 1-4 1074? Vem"!. 
The total potential difference across the current system is 1-4 x 104 V, a value large 
enough to be significant in auroral theory. 
The Pedersen current J, flowing during a moderate stor will be about half 
Jo or 10~4e.m.u. per cm strip. The total current flowing out of the cylinder is 
2aRJ, or 6-3 x 10*e.m.u.; this is also the total current flowing down the cylinder. 
In the case of the geomagnetic field the radius R of the bundle of lines of force 
varies with the distance from the Earth but the product 27R/J, is constant so that 
equation (21) gives 
ho = 1°3 x 10!R-?, (22) 


The angle through which a line of force rotates in a distance dz is 


ho 1*3 x 104 
dé =-— dz = — —dz. (23) 
H; R?H, 


The magnetic flux in the bundle is 7R2H, and must be a constant so that equa- 
tion (23) may be integrated to give 


13 x 1042 sia 
= —__—_., 2 
RH, . 


The value of R?H; at the Earth is about 5 x 1015 and if z is taken as say 6 earth 
radii then @ ~ 0-1 radian which provides a rough estimate of the total twist. If 
this is related to the form of distortion depicted in Figure 5, then the outermost 
lines of force will twist through ~ 0-2 radian and those much closer to the Earth 
through zero angle. ‘The twist relative to an intermediate region is then ~0-1 
radian. 

A twist of 0-2 radian originating in the region where hy, ~ Ho ~ 1073 gauss, 
is easily accounted for by a modest solar wind provided B in Equation (7) is not 
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too small. For example, if 8 = 0-2 and u = 108cms~! then the number density 
required is N ~ 50cm~% and the particle flux Nu is about 5 x 109 protons cm~? s~}. 
This is about 20 times less than some estimates of the peak value. 

The greatest uncertainty in the above estimate of the solar wind flux is the value 
of the “frictional” term 8. An order-of-magnitude estimate of B may be made by 
comparing the magnitudes of observed Dst and Ds. The former results from an 
isotropically propagating wave, the latter is focused by converging lines of force 
and so should increase in strength by a factor (H,/hm)', where H, is the strength 
of the main field near the Earth (~0-4 gauss). Putting hp, = 107% gauss as above 
the gain factor is 20 whereas the observed difference in intensity is about 4. This 
suggests that 8 might have a value of 0-4. There is another factor which may in- 
crease the transmission efficiency of the Ds wave, however, that is the fact that 
the observed disturbance is due to Hall current and stronger than the original 
current by a factor fogds + foids ~ 2 (Baker & Martyn 1953). This would reduce 
B to the value o-2. 


9g. Some secondary storm effects (including SC*) 

Of the various features of a Ds disturbance listed in Section 2 above, three 
remain to be explained, perhaps the most notable being the preliminary reverse 
impulse. 


g.1. The preliminary reverse impulse Ds (PR or SC*).—In Section 7 it was 
shown that if the current J; at the leading edge of the twist in a bundle of lines of 
force could not escape freely, then its reverse space-charge field would cause 


slipping. This effect is simply relative rotation of plasma and field as in an un- 
loaded homopolar generator. ‘The leading edge of the slip region moves with the 
Alfvén velocity but the arrival of the full twist and its radical current system J, 
is delayed. The arrival of the reverse electric field is not delayed however ; in fact 
it might even precede the wave front because of the flow of space charge down 
the lines of force. 

In the example given above and in Figure 6, negative space-charge develops 
inside the cylinder of lines of force and positive charge on its surface. The electric 
field is directed inwards and the ionospheric Hall current flows in an anticlockwise 
direction, opposite to that during the subsequent phases. In this way the two 
features of Ds(PR) are explained: first it is almost an exactly reversed (and 
weakened) Ds(SC). Second, the presence of Ds(PR) implies a delay in the onset 
of Ds(.SC) (Gerard 1959). 


9.2. Asymmetry about the noon meridian.—From the above theory, as devel- 
oped to this point, perfect symmetry of the Ds field about the noon meridian would 
be anticipated. However, instead of the centres of the equivalent current systems 
(A and B of Figure 1b) being at 06 hours and 18 hours local time they are an hour 
or so earlier. This time difference varies from storm to storm and by smaller 
amounts during a particular storm. 

The obvious asymmetrical feature in the interaction of the solar wind and geo- 
magnetic field is that the heavier protons are deflected to the right and the lighter 
electrons to the left. This feature might account for the positions of the Ds 
current loops in the following manner. In Section 5 it was stated that within the 
region of interaction of the geomagnetic field and the solar wind Hall current was 
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not of primary importance because the electric fields v AH and Ep, nearly balanced 
and so the ions followed nearly straight lines. However, they will not follow 
perfectly straight lines and any residual deflecting force will affect the electrons 
much more than the ions. 

If the electrons suffer appreciable deflection they will lag behind the ions, 
thus creating an electric current in the direction v of the solar wind. This current 
is simply Hall current and flows in both layers of gas outside the dashed lines (RR) 
of Figure 4. Since this current flows only in the gas passing through the geo- 
magnetic field, it is not continuous and causes positive space charge to accumu- 
late in the dark hemisphere (at PP) and negative charge in the light hemisphere. 
These accumulations shift the “‘centres of gravity” of the space charge on the 
boundaries RR in a clockwise direction both in the AM and PM hemispheres. 
This movement in turn moves the centres of the hydromagnetic twists and so the 
ionospheric current loops to positions of earlier local time. 


9.3. The low-latitude current loops of Ds; MP).—The remaining feature of Ds 
requiring explanation is the pair of weak current loops of Figure rb centred at 
latitude ~ 58°. ‘These loops and the squeezing of the main loop systems A and B 
up into the region in and above the auroral zone are the main differences between 
the Ds(SC) and Ds(MP) phases. 

These differences should be accounted for in terms of the different origins of 
Ds(SC) and Ds(MP) described above. The former is caused by the solar gas 
impinging on the solar side of the geomagnetic field and bending the lines of force 
out of their meridian planes. This bending away from the Sun causes a clock- 
wise ionospheric current loop in the AM hemisphere. As the main phase deve- 
lops, lines of force from all longitudes are swept away to form a geomagnetic tail. 
Only lines of force which cut the Earth at high latitudes contribute to this tail and 
after the solar gas cloud has passed the main twists and ionospheric current loops 
are associated only with these lines. Hence the current loops should be concen- 
trated at high latitudes as observed. 

Apart from the tail section, the geomagnetic field is now free from forces in 
the Sun—Earth direction. It will, however, experience forces from the opposite 
direction due to expansion of the bundle of lines of force forming the tail. For the 
present discussion this bundle may be regarded as consisting of three parts. The 
first is a central part whose lines of force cut the Earth at longitudes a few hours 
from local midnight. The pressure of this bundle in tending to expand has much 
the same effect as the pressure of the solar wind near the midday meridian. It 
causes a small increase in the horizontal component which partially neutralizes 
the general decrease (Dst main phase) resulting from the formation of the whole 
tail (P.2). 

The two outer parts of the geomagnetic field comprise lines of force which have 
been pulled away from the sunlit hemisphere. In tending to expand they push 
those remaining parts of the geomagnetic field lying near the 06 and 18 hour 
meridians in the direction of the Sun. This pressure should have the reverse effect 
of that of the solar wind blowing past these same lines of force (Ds(.SC)) and so 
should lead to reverse current loops at latitudes below the auroral zones. 


10. Discussion 


In Section 3 and Appendix it was shown probable that the equivalent current 
systems of all three phases of a polar geomagnetic storm are real current systems 
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flowing in the ionosphere and that the main observed effects are due to Hall 
currents. These quasi-steady currents may only be maintained by electric cur- 
rents flowing into and out of the ionosphere from the region of interaction of the 
solar wind and the geomagnetic field. 

The origin and propagation of those external disturbances have been studied 
quantitatively and found to fit the observational results. Essentially they are due 
to a type of hydromagnetic disturbance which may only travel along the magnetic 
lines of force. The important effect during the main phase of a storm is the forma- 
tion of a geomagnetic “tail” directed away from the Sun. 

Although the discussion has centred on sudden commencement storms because 
these comprise such a variety of features, it applies equally well to the storms 
which start slowly, last for 3-4 days or more and have a 27-day recurrence period. 
The explanation of these storms is similar to the Ds(SC) or Ds(MP) above. 

It is important to notice that although storm effects have been arbitrarily 
divided into longitude independent (Dst) and longitude dependent (Ds) parts any 
satisfactory theory of the main phase must explain both parts. It would seem that 
any ring-current theory or any theory depending on trapped (van Allen) radiation 
could only account for circularly symmetrical effects (Dst) while the magnetic tail 
theory introduces the important element of asymmetry. 


Appendix 


There is now fairly conclusive evidence that the various geomagnetic storm 
effects result, indirectly, from temporary increases in the flow of solar ions past the 
geomagnetic field. If the direct cause of the storms (or of their polar parts) is a 


system of ionospheric winds then the latter must be caused by the solar wind. This 
could only be effected by hydromagnetic action and if this occurs we already have 
a source of magnetic disturbance without postulating violent ionospheric wind 
systems. To be more specific it has been shown that geomagnetic disturbances in 
general pass through the ionosphere as through a solid semi-conductor (P.1). 
This is because neutral atoms collide so infrequently with ions that it takes a 
number of hours (~6) to start ionospheric winds by means of electromagnetic 
forces. It is concluded that a storm theory based mainly on changes in ionospheric 
winds is untenable. 

Another form of atmospheric dynamo theory is that of Obayashi & Jacobs 
(1957) in which the wind system remains that of geomagnetically quiet days but 
the electrical conductivity (and hence the current density) increases. This theory 
and the earlier version (or their combination) has been criticized by Maeda (1959) 
on the grounds that the winds or conductivity changes needed would be un- 
reasonably large. The quiet day ionospheric current system driven by the wind 
dynamo flows in low as well as high latitudes and it is known that no substantial 
increase in the conductivity at these low latitudes occurs at the time of a storm. 
An examination of the whole dynamo current system suggests that even if the con- 
ductivity near the poles became infinite the current around the whole dynamo 
circuit would show no great increase (D. F. Martyn, private communication). 
Added to these difficulties is that of the Ds(PR) which shows a completely reversed 
system of currents and so could hardly be accounted for in terms of changing 
conductivity. 


It is concluded that ionospheric winds play no important part in driving the 
currents responsible for Ds. 
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There is another way in which the electric field E; of Figure 2 could be induced, 
apart from ionospheric winds flowing radially inwards towards O. That is, by an 
outward movement of magnetic lines of force and hence a decrease in H (equa- 
tion (1)). Such movement would require a radially outward stress on the bundle 
of lines of force above the current loop. In view of the results of P.1 on hydro- 
magnetic transmission in this region it seems most unlikely that the plasma there 
could continue to exert such a mechanical force for the many hours of Ds(SC) 
and Ds(MP). However, even if it could, it is easily shown that the low conduc- 
tivity of the ionosphere would allow such a large value of —dH/ét that the 
vertical component would decrease much more than observed. This would mean 
that the inductive current loop opposing the decrease in H (which is mainly vertical 
component) would transfer to the Earth itself. Here the conductivity is sufficient 
but the components of the perturbation field other than the vertical component 
are quite different from those inferred from Figure 1(b). 

Thus the electric field driving the Ds current system cannot be an induction 
field (irrespective of where the currents flow) and the current system cannot flow 
in the Earth, where only an induction field would be effective. 


C.S.I.R.O., 
Radiophysics Division, 
Sydney, 
Australia: 
1960 April. 
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Micropulsations 
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Summary 


A detailed investigation is carried out on geomagnetic micropulsa- 
tions recorded by rapid-run magnetograms during the International 
Geophysical Year from 17 observatories as widely distributed geo- 
graphically as possible. A number of well-defined cases (14 Pt’s and 
3 Pc’s) are investigated in detail. Diurnal and world-wide charac- 
teristics are derived and the equivalent overhead current systems 
which may give rise to the micropulsations are constructed. 

Many of the examples of Pt’s can be divided into two groups: viz. 
Pt- which is a Pt associated with a negative bay in the auroral zone 
and Pt+ which is a Pt either associated with a positive bay or one 
which appears with no bay in the auroral zone. Although both groups 
of Pt’s, with their associated bays, have different characteristics, there 
are many points of similarity. Both classes have more or less coherent, 
synchronous wave forms over a wide area, i.e. they have the same period 
at different stations, where they are either nearly in or out of phase. 

It is found that Pc’s consist of two different wave bands, both 
occurring during the daylight hours. The first band contains shorter 
period Pc’s (15-308) with maximum amplitudes of a few gammas in 
moderate latitudes. They are not synchronous over wide areas. The 
second band contains longer period Pc’s (30-gos) with maximum 
amplitudes in the polar regions, amounting, on occasions, to as much 
as a few tens of gammas. They are well synchronized round the 
world. 

Attention is drawn to the occurrence of long period, continuous, 
pulsations (called LPc’s) which have a longer period than Pec’s and 
usually appear simultaneously with Pc’s in polar regions. 

The characteristics of micropulsations found in this investigation 
may throw some light on their origin. In this regard, a tentative dis- 
cussion is given of hydromagnetic oscillations in the outer atmosphere 
which are excited by solar corpuscular streams. 


1. Introduction 


Geomagnetic micropulsations are fluctuations of the Earth’s magnetic field 
whose amplitudes range from a fraction of a gamma to, on rare occasions, as much 
as a few tens of gammas. Periods range from about 0-18 to 10min. Most work 


* On leave of absence from the Radio Research Laboratories, Kokubunji, Tokyo, Japan. 
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has been done in the lower frequency end of this range, very little being known 
at the higher frequencies from o-1 to 1oc/s. In recent years, mainly as a result 
of the greatly increased number of continuous observations that are being made, 
certain definite types of micropulsations have been recognized and some systematic 
classification has been attempted. Ironically enough the great increase in the 
amount of data obtained during the International Geophysical Year has only 
served to show in some cases that the earlier classification was too simple and that 
in reality there exist many more different types or variations of the earlier, broader 
classification. ‘This classification divided micropulsations into three categories, 
Pc’s, Pt’s and Pg’s, which are defined below. 


Pc’s (pulsations continues) 


A series of pulsations lasting for many hours with periods usually in the range 
of from 10 to 60s and amplitudes of the order of o-1 y. The maximum occurrence 
frequency is during the morning hours. Well-defined oscillations with periods 
less than 10s have been observed by Troitskaya (1959a), Benioff (1960) and other 
workers. There is no doubt that they form a distinct class from Pc’s, but they are 
not investigated in this report. 


Pt’s (trains de pulsations) 


These appear as several series of oscillations, each series usually lasting from 
10 to 20min the whole phenomenon lasting for periods of not more than about 
th. They are well damped with longer period than Pc’s (from 40s to a few min) 
and usually of greater amplitude, i.e. of the order 0-5 y. The maximum occur- 
rence frequency is before midnight. 


Pg’s (giant pulsations) 


A series of pulsations of large amplitude (up to a few tens of gammas) appear- 
ing only in or near the auroral zones. The period is longer than that of Pc’s (up 
to several minutes), and the duration is of the order of an hour or more. 

The results of several Soviet investigators were presented at a recent sym- 
posium on rapid variations in the Earth’s field held under the auspices of the 
International Association of Geomagnetism and Aeronomy at Utrecht, 1959 
September. Mainly as a result of their work in the Arctic and the Antarctic, 
Troitskaya (1959b) recommended that Pc’s be subdivided into three groups, with 
periods in the range 5~15s, 20-40s and 50-gos. Of particular interest is the sharp 
fall in the number of hours that recorded Pc’s (in the 20-30s range period) 
during the middle months of the long polar night. There are also differences in 
the characteristics of Pt’s in polar regions and middle latitudes. ‘Troitskaya 
stressed the importance of shorter period pulsations, which fall into two 
classes—Pc’s with periods from 5—15s, characteristic mainly of the polar regions 
and whose diurnal distribution is quite different from that of Pc’s of longer 
period (20-308), and beating pulsations (pearls) with periods from 1-48. 

Kalashnikov (1959) maintained that the present classification Pc and Pt is too 
subjective and that definitions of pulsations should be based on more objective 
elements, period, amplitude and duration of the oscillations. Series of oscillations 
of different amplitude and period are often superimposed and may be caused by 
different sources and different excitation processes. Kalashnikov stated that at 
the Soviet stations there were often periods when the field of rapid variations 
became non-homogeneous and noticeable differences of the order of tenths of a 
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gamma existed at a distance of 1 km both in latitudinal and longitudinal directions. 

It was to investigate further some of these problems that a detailed analysis 
was carried out of a number of well defined cases of Pc’s and Pt’s, using data 
from stations as widely distributed as possible. This enabled the equivalent over- 
head current systems, which may give rise to these micro-pulsations, to be 
constructed. 


2. Source of data 


In order to investigate the morphology of micropulsations, it is desirable not 
only to obtain information on their average characteristics, but also some know- 
ledge of their world-wide behaviour. Unfortunately, magnetic stations, operating 
rapid-run magnetograms, are not well distributed geographically. The only reason- 
able coverage is in North America and the Pacific, the stations being mainly 
operated by the United States Coast and Geodetic Survey. The observatories 
whose data are used in the present study are given in Table 1. 


Table 1 
Stations supplying data used in this report 


Geographic Geomagnetic 


~ 


Station Abbrevia- Lat. Long. Lat. Long. 
tion 

Point Barrow PB 71° 18’N 156° 46’W 68°-6N 241°-0E 
College Co 64° 51’N 147° 50°W 64°-7N 256°-sE 
Big Delta BD 64° o9’N 145° 51°'W 64°-4N 259°°0E 
Healy 63° 51’N 149° 00’ W 63°-6N 256°°sE 
Sitka Si °03’N 135° 20°W 60°-oN 275°'4E 
Lovo °2a1’N 17° 50°E 58°-2N 105 8E 
Fredericksburg , > 12’N 77° 22’'W 49°°6N 349°°9E 
Tucson > 114’N 110° 57°W 40°°4N 312 
Kakioka < ° 14’N 140° 11’E 26°-oN 206°: 
Honolulu 18’N 158° 08’ W 21°-oN 266°: 
Guam ; °27’N 144° 45'E 3°-9N 212°'8 
Koror > 16'N 134° 32'E 3°°38 203°°5 
Apia J > 48'S 171° 46°W 16°°0S 260°: 
Hermanus > 24'S 19° 13'E 33°°3S 80°; 
Byrd Station y 59'S 120° o1'W 70° 6S 336°° 
Little America LA 18'S 162° 10°W 74°°0S 312°-0E 
Wilkes Wk * 15°85 110° 31’E 77°°8S 179° 0E 


The data from these 17 observatories were recorded photographically, not 
only on rapid-run magnetograms with a recording speed of 20mm per 5 min, but 
also on standard speed magnetograms with a recording speed of 2omm/h. Addi- 
tional data were obtained from the publications of the International Association of 
Geomagnetism and Aeronomy (supplied by Dr J. Veldkamp, Kon. Nederlansch. 
Meteorologisch Instituut, De Bilt, Holland). 

A number of well-defined micropulsations, which occurred mainly in American 
regions (3 cases of Pc’s, 14 cases of Pt’s and 5 cases unclassified) are considered. 
After reading and comparing the data, diurnal and world-wide characteristics are 
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derived, and the equivalent overhead current systems drawn for some typical 
cases. ‘Table 2 gives the times and dates of the examples of micropulsations which 
are analysed in this paper. The table also contains information on the newly 
defined pulsation, LPc, which appears mainly in polar regions, with long-period 


(3-7 m) continuous oscillations, lasting several hours. 


3. Characteristics of Pc’s 


(a) Amplitude distribution.—Pc’s are defined as continuous pulsations usually 
appearing on the daylight side of the Earth, lasting for many hours with periods in 
the range from 10 to 60s, and amplitudes of the order of a few tenths of a gamma 
in middle latitudes. However, if the evidence of all investigators over a wide area 
is considered, it becomes apparent that Pc’s do not correspond toa single wave band 
but to two wave bands, the first a shorter period Pc (about 15—30s) and the second 
a longer period Pc (about 30-gos). Figure 1 shows the spectrum of Pe’s which 
have been investigated by several workers from stations all over the world, and 
clearly illustrates the division into these two different wave bands. 


}e—1AGA Resolution of period of Pc’s ~ 


Shorter Pc Longer Pc 


- 


— ce 60 io 
Period of Pcs (Seconds) —— 
Fic. 1.—Latitudinal distribution of spectrum of Pc’s with periods 
greater than 10s as reported by different investigators from data shown 
on the right-hand side. 


In a previous investigation (Jacobs & Sinno 1959, 1960), latitudinal and 
longitudinal characteristics were obtained of the diurnal occurrence frequency of 
Pc’s during the International Geophysical Year. It was found that the occurrence 
of Pc’s increased from the equator towards higher latitudes regardless of period. 
The present investigation shows that for one type of Pe (the longer-period Pc), the 
amplitude also increases from the equator towards the poles. This is shown in the 
examples observed on 1958 February 12 (Figures 2 and 3) and 1958 April 2 
(Figure 2), which have periods of 33s and 45s, respectively. These examples of 
Pc’s have coherent, synchronous wave forms over a wide area. As can be seen from 
these figures, Pc’s of this type often have amplitudes up to ten gammas or more 
in polar regions. 

A pulsation appeared on 1958 May 7 (Figures 2 and 4), with a period of about 
20s and with maximum amplitude at about latitude 50°. Many similar examples 
(such as those which occurred on 1958 February 8 and May 28) have been found: 
they all have in common comparatively short periods (from about 15 to 30s). 
These shorter-period Pc’s show different periods at different latitudes although it 
is not possible to be definite on this point since the speed of the rapid run record- 
ings is too slow for any fine resolution. Although wave bands below 15s period 
have been defined by several investigators (Troitskaya, Maple, Benioff, Duffus, 
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etc.) they have not been included in the present report because the records were 
obtained at too slow a speed to identify them. 

According to Maple (1959), oscillations with a period of about 20s occur more 
often than those of longer period (70s) at Tucson (geomagnetic latitude 40° 4’ N). 
On the other hand, 70s oscillations predominate at Eskdalemuir (geomagnetic 
latitude 58° 4’ N). This frequency dependence ties in with the results of the present 
investigation, where short-period Pc’s were found to have their maximum amplitude 
at comparatively low latitudes. The latitudinal distribution of amplitude for the 
short-period Pc’s is quite different from that of the longer-period Pc’s which have 
maximum amplitudes in the polar regions, the amplitudes decreasing at lower 





1957 Sept. 23 


1958 Apr. 2 





— . 
72 lUlfOe hUO 5O O 
Geomagnetic Latitude 











Fic. 2.—Amplitude distributions for various Pc’s which occurred on 
1957 September 23, 1958 February 12, April 2 and May 7. 


latitudes. This tendency is also shown in Figure 1 where the period at the lower 
boundary of the longer Pc wave band increases with decreasing latitude. More- 
over, of the two examples of longer-period Pc’s, that with the shorter period has a 
greater range of amplitude variation as is shown in Figure 3. Combining the 
results of previous investigations with the latitudinal distribution of the amplitude 
of Pc’s which has been derived in the present analysis, equi-amplitude contours 
can be drawn as shown in Figure 5. 

Other micropulsations, which could not be rigorously classified or analysed, 
were observed. Thus a pulsation appeared on 1957 September 23, with a very 
enhanced amplitude—up to several tens of gammas in the auroral zone (see 
Figure 2). However, this pulsation may not belong to the category of Pc’s, because 
of its irregular wave form and different period at distant observatories. It may form 
part of a series of irregular variations during a severe magnetic storm. It is worth 
noting, however, that the diurnal (longitudinal) characteristics of all types of Pc 
show a more or less daytime maximum as shown in Figure 6. 

Finally, attention is called to beating Pc’s which have been reported by Kato & 
Watanabe (1957, 1958). These are of longer period than the “pearls” of Troitskaya 
(1953, 1959a) and Benioff (1960). Beating Pc’s usually show a pulsation whose 
amplitude varies with a period of several minutes. Sometimes the modulating 
period is identical with that of LPc’s which often appear in polar regions. 

Y 
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(b) Equivalent overhead current system.—As shown in Figure 3, cases can be 
found where Pc’s occur simultaneously over a wide area. Strictly speaking, the 
recording paper speed of the rapid-run magnetogram is too slow to resolve the 
phase difference between distant stations. However, it is possible to determine 
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Fic. 3.—Traces from rapid-run magnetograms of longer-period Pc’s at 
various stations on 1958 February 12. 


whether the records are nearly in or out of phase. With this proviso, the ideal 
overhead current system which would induce the observed Pc, can be drawn as 
shown in Figure 7. The centre of the current vortex is located at about 50 degrees 
latitude and at about the sub-solar point. The concentrated current lines over the 


polar regions do not close in one sheet; these open-ended current contours may 
be extended to three-dimensional space. 
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Fic. 4.—Traces from rapid-run magnetograms of shorter-period Pc’s 
at various stations on 1958 May 7. 
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Fic. 5.—Equi-amplitude contours of Pc’s, derived from Figure 1. 
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Fic. 6.—Local time (longitudinal) variation of various Pc’s. 














Fic. 7.—Equivalent overhead current system for Pc at 2220-2225 U.T. 
1958 April 2. Electric current between successive lines is 5 x 10° A. 
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The equivalent overhead current system cannot be drawn for the shorter- 
period Pc’s because they are not synchronous at different stations. These shorter- 
period Pc’s may be caused by hydromagnetic waves in the lower atmosphere or 


ionosphere where local inhomogeneities, not present in the outer atmosphere, may 
develop. 


4- Characteristics of Pt’s 


(a) Amplitude distribution.—Pt’s are defined as well-damped transient oscilla- 
tions with periods longer than Pc’s. Generally Pt’s occur during the night time 
and often accompany a geomagnetic bay. In the present investigation the lati- 
tudinal dependence of the amplitude of several well-defined Pt’s is obtained. 
Fourteen examples of Pt’s with large amplitudes (see Table 2) were selected from 
data obtained from United States geomagnetic stations. 


Table 2 


Times of the occurrence of micropulsations analysed in this report 


Commence- Associa- 
Date ment Duration Type tion 


(U.T.) 


1957 Sept. 0402 6m Pt 

1957 Sept. 0800 4°5h LPc* 

1957 Sept. 0700 unclassified 

1957 Sept. 0440 Pt 

1957 Sept. 0456 Pt 

1957 Nov. 1555 LPc* 

1957 Dec. 0332 Pt 

1957 Dec. 0448 Pt 

1957 Dec. 25 0704 Pt 

1957 Dec. 0936 Pt severe bay 
1958 Jan. 0849 Pt severe bay 
1958 Feb. OSI c Pt 

1958 Feb. 0455 Pt 

1958 Feb. 2100 LPc* Pc 
1958 Feb. 2250 Pc LPec* 
1958 Feb. 1730 . LPc* 

1958 April 3 1800 Pc LPc* 
1958 May 1400 3 Px 

1958 May 0505 8m Pt severe bay 
1958 May 2339 Pt 

1958 June 0340 30m Pt 

1958 July 3 0640 35m Pt 


* Defined in this report. 


It was found that the latitudinal distribution of the amplitudes of Pt’s is 
different between that group of Pt’s associated with severe negative bays in the 
auroral zone and that group of P?’s associated with positive bays or without any 
accompanying bay. These two groups of P?’s are distinguished, as Pt- and 
Pt+, respectively. The latitudinal amplitude distributions derived from two 


examples are shown in Figures 8 and 9. The main features of these graphs are 
confirmed in many other cases. 








342 J. A. Jacobs and K. Sinno 


Pt~’s show a very sharp maximum in the auroral zones where the peak value 
reaches about ten gammas or more. On the other hand Pt+’s show a comparatively 
flat distribution with a maximum amplitude in the sub-auroral zone and subsidiary 
peaks in the polar and equatorial regions. These latitude trends of the amplitudes 
of Pt~’s and Pt*’s are the same as those of their associated bays, and lead to the 
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Fic. 8.—Amplitude distribution of Pt- and its associated bay, com- 
mencing 0936 U.T. 1957 December 26. 
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Fic. 9.—Amplitude distribution of Pt*+ and its associated bay com- 
mencing 0051 U.T. 1958 February 9. 


conclusion that Pt’s are fundamentally very closely related with bays. Thus the 
causal mechanism of these two phenomena may well be identical. 


(b) Equivalent overhead current system.—The equivalent overhead current 
systems for an example of Pt- and Pt* are shown in Figures 10 and 11 respec- 
tively. Although both types of Pt’s have synchronous wave forms at different 
stations, it is not easy to draw the contours in the case of Pt- on account of the 
complexity of the wave form, especially near the auroral zones. For simplicity, 
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the direction of the current contours and magnetic arrows are shown for the maxi- 
mum amplitude of the oscillation which fluctuates with the period of the P?’s. 
This procedure was also adopted in drawing the overhead current systems for Pc’s. 

Both the contours for Pt- and Pt+ resemble the well-known bay current 
systems which have been drawn by Silsbee & Vestine (1942), although Pt’s have 
only one intense current vortex centre, on the night side. The current vortex 
centres are located at about 50—60 degrees latitude which is in the sub-auroral zone. 














Fic. 10.—Equivalent overhead current system for Pt- at 0936 U.T. 
1957 December 26. Electric current between successive lines is 
5 x 107 A. 


5. Long-period continuous pulsations (/Pc) in polar regions 


(a) Amplitude distribution.—Figure 12 shows part of the traces of standard 
magnetograms, speed 20 mm/h, from College, Alaska, and Fredericksburg, Virginia, 
on 1958 February 12-13. Between 2300-2400 U.T. on February 12 series of very 
intense longer-period Pc’s were observed. They continued for about one hour 
and appear as thick traces on the magnetograms because of the slow paper speed. 
The same series of Pc’s were shown in Figure 3 recorded on rapid-run magneto- 
grams. Of special interest are the pre- and post-oscillations of this particular Pc. 
Figure 12 shows clearly the intense pulsation of about six minutes period which 
prevailed during the previous two hours and after the Pc. Usually this type of 
pulsation is observed in polar regions and is frequently associated with Pc’s; it has 
been called a long-period continuous pulsation (LPc) in this paper. 
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An example of an LPc on rapid-run magnetograms (which is also seen in 
Figure 3 as a pre-Pc pulsation), is shown in Figure 13. The oscillations, of about 
six minutes period, are nearly synchronous over a very wide area, from the Arctic 
to the Antarctic and in equatorial regions. Moreover, the polarization is nearly 
linear at each station. Sometimes, however, elliptically polarized LPc’s are 
observed as shown in Figure 14. The latitudinal dependence of the amplitude 
distributions for both Pe’s and LPc’s are very similar, as is shown in Figure 15. 
Usually the amplitudes of LPc’s are about ten times greater than Pe’s in polar 
regions. 














Fic. 11.—Equivalent overhead current system for Pt*+ at oo51 U.T. 1958 
February 9. Electric current between successive lines is 5 x 103 A. 


(b) Equivalent overhead current system.—Although the polarity of the horizontal 
component is not always exactly linear, small phase differences have been neglected 
for simplicity. The polarization when non-linear appears as a counter-clockwise, 
flat-shaped ellipse (see Figure 14). The equivalent overhead current system for 
LPc’s (Figure 16) is very similar to that for Pc’s except over the polar cap where 
LPc’s are not so well defined, as is shown, for example, in the traces at Point 
Barrow (see Figure 13). 


(c) Relation between long-period Pc’s and LPc’s.—Figure 13 shows the remark- 
able time association of Pc’s and LPc’s, although LPc’s do not usually occur in 
middle latitudes. This is further illustrated in Figure 17 which shows the occur- 
rence of Pc’s and LPc’s at Point Barrow, Alaska, during 1958 February. The 
diurnal occurrence of LPc’s, like that of Pc’s, has a daytime maximum a few hours 
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Fic. 12.—Examples of Pe and LPc. Traces are from ordinary magneto- 
grams at College and Fredericksburg on 1958 February 12-13. 
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Fic. 13.—Example of LPc. Traces from rapid-run magnetograms at 
various stations on 1958 February 12. 
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before local noon. Besides this time association there is further evidence of the 
relation between Pc’s and LPc’s. Sometimes beating oscillations of Pc’s, which 
appear mainly in middle latitudes, occur simultaneously with LPc’s in polar regions. 
Moreover, the periods of both the beating Pc’s and the LPc’s are about the same, 
viz. from 3 to 7m (Kato & Watanabe 1957, 1958). These close relations between 
Pe’s and LPc’s indicate that both pulsations have a common source, probably outer 
atmospheric oscillations. 


6. Discussion of the results 


(a) Hydromagnetic oscillations in the outer atmosphere.—Although equations for 
hydromagnetic oscillations of the outer atmosphere have been set up by Dungey 
(1954), they are too complex to be solved in the general case. However, if there is 
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Fic. 14.—Example of LPc showing the polarization of the horizontal 
component. ‘Traces from rapid-run magnetogram at Sitka on 1957 
November 24. 


no coupling between poloidal and toroidal oscillations, solutions may be obtained 
from the following two equations: 
Poloidal 

aires hdl “(sin dy a i — — —\(Rsind Rng @ 
oH? ot <i 


Toroidal 


p & ; uy 
P — _ (Rsin6)-2[(Ho . ; 
Ee ot? _— eal : (2) 


Spherical polar coordinate (R, 0, ¢) are used where R is the distance from the 
centre of the Earth, p the plasma density of the outer atmosphere, and po the 
permeability (mks units are used). The magnetic field H is considered to be the 
sum of a constant dipole field Hp and a small disturbance h introduced through 
the motion of the medium, where |Ho| > |h]. 

The solutions of equations (1) and (2) involve the following components of the 
electric field E, velocity u and disturbed magnetic field h. 


Poloidal 
E(o, Oo, E,), u(up, Up, o), h(hpr, ho, 0) (3) 

Toroidal 
E(En, Ey, 0), w(0, 044), (0, 0, hy). (4) 


The toroidal mode was studied by Dungey (1954), Kato & Watanabe (1956), 
and Obayashi & Jacobs (1957) as the possible origin of giant pulsations (Pg’s), 
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since the toroidal mode and Pg’s show similar characteristics. In particular the 
occurrence of Pg’s is restricted to places near the auroral zones and the main oscil- 
lating component is the D-component of the magnetic field. On the other hand, 
the synchronous nature over a wide area and the approximate meridional plane 
of polarization of Pc’s and Pt’s might be explained by the poloidal mode, as has 
been pointed out by Kato & Akasofu (1956) and Kato & Watanabe (1957, 1958). 
Following their analysis, and assuming that the field quantities depend on t through 
a factor exp(iwt), equation (1) can be written in the form 


/ gz 1—u2) 22 “ 
(—+! = a) - (5) 
where 

Rsind.E, =, cosé =p (6) 
and 


pw?/uoHo® = 1/7). (7) 


Writing ¢ = fi(x) . fo(x) where R = yx, equation (5) can be solved by the method 
of separation of variables, f; and fe satisfying the equations 


os (1- n(n +1) 


72 )f =o (8) 


x2 


d*f. 
G-n pt n(nt 1)f2 = 0. (9) 
be 


For any integral value of n, the exact mode of the oscillations may be obtained 
from equations (8) and (9), if the variation of the Alfvén velocity V4 = Ho(yo/p)+ 
= wy in the outer atmosphere is known. However, without knowing the exact 
boundary conditions, it is possible to make a rough estimate of the period of 
the standing wave in the outer atmosphere which is confined between its inner 
and outer boundaries. 


From the simple physical analogy of standing waves, the fundamental period 
and its harmonics are given by 


(10) 


where 7, is the period of the harmonic of degree n and 7; the fundamental period. 
r; and re are the radii of the inner and outer boundaries of the outer atmosphere 
measured in units of the Earth’s radius a. This relation implies eigen periods of 
the outer atmospheric oscillation equal to twice the travel time between the outer 
and inner boundaries. On the basis of our present knowledge, 7} is about a few 
hundred seconds. This period is too large for Pc’s and Pt’s to be the direct expres- 
sion of the fundamental wave. 

Recently Dessler (1958, 1959) has shown that there are sharp gradients of the 
Alfvén velocity at the top and bottom of the outer atmosphere (as shown in the 
right half of Figure 18) where the waves are partially reflected back. Besides these 
sharp boundaries, equation (8) shows the existence of another barrier at the height 
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Fic. 15.—Amplitude distributions of Pc and LPc at 2215-2230 U.T. 
1958 February 12. 

















Fic. 16.—Equivalent overhead current system of LPc at 2215-2230 U.T. 
1958 February 12. Electric current between successive lines is 
5 x10* A. 
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Fic. 17.—Hourly occurrence of Pc’s and LPc’s at Point Barrow in 
1958 February. 
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of the Alfvén velocity maximum which becomes effective when the period of the 
waves is longer than about 10s although it may not become complete until the 
period exceeds 100s. This barrier for an Alfvén wave is similar to that caused by 
the ionosphere for an electromagnetic wave. 

Although the exact behaviour of the propagation of Alfvén waves in the outer 
atmosphere cannot be established without further detailed knowledge of the outer 


atmosphere, plausible mechanisms for the generation of Pc’s, LPc’s and Pt’s are 
given below. 


(b) Model for Pc’s and LPc’s.—In the outer atmospheric cavity standing waves 
whose periods are longer than about a minute cannot leak downwards sufficiently 
to be observed because the barrier is so high and thick. Thus only the higher 
harmonic waves 7, whose periods are less than a minute would be observed. 
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Fic. 18.—Schematic representation of Pc waves trapped in the outer 

atmosphere. The boundaries of this standing wave are formed by the 

sharp gradient of the Alfvén wave velocity shown on the right-hand side 
of the figure. 


These higher harmonic waves interfere with each other and set up beating oscilla- 
tions which are modulated by the fundamental period (Kato & Watanabe 1957, 
1958). The beating of the higher harmonic waves is represented schematically on 
the left-hand side of Figure 18. The evidence that Pc’s sometimes show a beating 
oscillation with a period of about five minutes can be interpreted by the above 
model. LPc’s which appear mainly in the polar regions might be an expression 
of the leakage of the fundamental mode under the special circumstances that exist 
in the polar outer atmosphere, such as the remarkable difference in temperature 
between the auroral zone and moderate latitudes, as has been observed by Horo- 
witz & LaGow (1958). The maximum occurrence of Pc’s around local noon can 
be understood since the perturbing force is located on the solar hemisphere where 
the corpuscular streams meet the outer atmospheric boundary. The centre of 
the current vortex in the equivalent overhead current system for the longer-period 
Pc’s at about latitude 50° and the latitudinal distribution of amplitude show that 
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the coefficients of the spherical harmonics predominate at nm = 3, 5 and 7. The 
odd numbers come from symmetry with respect to the equator. 

Since it is non-synchronous, it is not possible to consider the shorter-period Pe 
as originating as an outer atmospheric poloidal oscillation. It may, however, arise 
from hydromagnetic oscillations of the lower outer atmosphere or ionosphere 
below the barrier, since it shows local characteristics. The positive correlation 
between the occurrence of the shorter-period Pe and geomagnetic activity, and 
the approximate negative correlation of the occurrence of the longer-period Pc 
(Maple 1959), can easily be understood. Geomagnetic activity may help to develop 
the barrier, which would promote oscillations below, and at the same time inhibit 
the wave meeting the barrier from above. 


(c) Model for Pt’s—The maximum occurrence of Pt’s at night and in the 
auroral zones leads to the conclusion that they are outer atmospheric oscillations 
which are excited by the same charged particles that cause aurora and bays. 
However, in spite of this close correlation, aurora and bays occur more frequently 
with increasing solar activity during the sun-spot cycle whereas the converse is 
true for Pt’s. This inverse relation between the occurrence of Pt’s and solar 
activity has been reported by Yanagihara (1959). It might be understood if the 
barrier were to develop to such an extent at the maximum of solar activity that 
the waves responsible for Pt’s were inhibited, in spite of the increased frequency of 
aurora and bays. It is reasonable to consider that the barrier may depend not only 
on geomagnetic activity (as is the case for Pc’s), but also on solar activity. 

Although no typical Pg’s were observed during the International Geophysical 
Year, the similar inverse relation between solar activity and Pg’s, which is well 
known (Sucksdorff 1939), may be explained by the same reasoning; the barrier 
will exist for toroidal as well as for poloidal oscillations. 
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Note on Cusps in Seismic Travel-times 


K. E. Bullen 
(Received 1960 June 20) 


Summary 


When there is triplication in travel-time curves, it is often assumed 
that the cusps in the curves are necessarily associated with large wave 
amplitudes. Models which give a good representation of various types 
of seismic velocity variation are here set up to show that there are 
important cases in which this assumption does not apply. The analysis, 
which is facilitated by use of the variable « = 2d log r/d log n, shows 
fairly directly that dA/dp, where A and p are the angular length and 
parameter of ray, is sometimes discontinuous, sometimes continuous, at 
a cusp. 


1. Preliminaries 


In this note, an examination of the character of cusps in seismic travel-time 
curves is carried to greater detail than previously (Bullen 1945, 1959). It has 
often, but erroneously, been taken for granted that a cusp in the graph of the 
travel-time ¢ against the angular distance A is necessarily associated with abnormally 
large wave amplitudes. A purpose of the present note is to help clear up the 
confusion which appears to exist on this question. 

Let wv be the (P or S) seismic velocity at distance r from the centre of an Earth 
model (as usual, taken to be spherically symmetrical). The question at issue can 
be sharply focused by representing relevant features of velocity distributions 
in terms of » and a, where 

7 = 1/0, (1) 
a = 2dlogr/dlog». (2) 


The utility of the variable « has been previously demonstrated (1945, 1960), 
especially in respect of the quality of the approximation given by taking « constant 
over ranges of values of r. 

Cusps may arise in t—A graphs when v or its gradient (—dv/dr) increases 
discontinuously with increase of depth, or when the gradient increases con- 
tinuously but rapidly. Each of these cases will be considered. 

By (1), a positive jump in v corresponds to a negative jump in 7. By (2), a 
positive jump in —dv/dr corresponds to a negative jump in «a. 

We shall consider families of rays terminating at both ends at a level-surface 
r = ro. Other level-surfaces, r = r; and r = re, where re < 7 < 70, will also be 
referred to. The subscripts 0, 1 and 2 will indicate values of functions of r at 
r = 70, 71, Y2, respectively. 
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2. First model 


We consider first a model containing an outer layer L, for which 7) < r < 19 
and in which « takes the constant value ao, overlying a layer N, for which r < 1, 
and in which « = cao, where « is constant. The velocity v is taken to be con- 
tinuous at the boundary r = 7;, while the velocity gradient is taken to jump dis- 
continuously from LZ to N; thus 7 is continuous at the boundary, and « < 1. 

For rays entirely in L, we have (1959, section 7.3) 


A = agcos—!(p/no), (3) 
t = a9(y0?—4°)', (4) 
dA/dp = — (yo? — p?)*, (5) 


where p is the usual ray-parameter. 
For rays penetrating N, it can be deduced (fuller detail is being published in a 
further paper on ray theory) that 


A = ap{cos“(p/no) — (1 — x) cos~1(p/m1)}, (6) 
t = ao{(n0?— p*)! —(1 — «)m?— pi}, (7) 
dA/dp = ao\ — (0? — p?)* + (1 — «)(m? — p?) 4}. (8) 


Consider variations as p decreases continuously from the value yo, for which 
A =o. At p = x, dA/dp is by (5) equal to — «0; dA/dp then increases continuously 
to the moderate value — ao(yo?—p”)-* at p = m, where, by (8), it jumps discon- 
tinuously to + 0 (since x < 1). As p further decreases, (8) shows that dA/dp 
decreases continuously from + 0, passing through a zero at p = 7’, where 


m?—1'2 = (1—«)*(q0?—72). (9) 


Inspection of (3), (4), (6) and (7) shows that each of A and t¢ is continuous at 
both 4 and 7’. 

From these details, it is easy to construct Figures 1, 2 and 3 in which the 
points A, A’, A’’ correspond to one another, etc. Figure 3 shows the well-known 
triplication of the time-curve for this class of case. 

Figures 1 and 2, however, reveal an important feature not made explicit in the 
writer’s previous work (1945). Whereas C’ corresponds to an ordinary minimum 
of A (against p), at which dA/dp is continuous, B’ corresponds to a node-point 
maximum at which dA/dp is discontinuous and does not in general vanish. ‘Thus 
there is a significant difference in the character of the cusps at B and C which is 
not discernible from Figure 3 alone. 

The difference is of practical importance, since p is proportional to cos e, 
where e is the angle of emergence of a ray, so that a condition for large wave 
amplitudes at a given value of A is dA/dp = o. This condition is realized at C 
but not at B. Thus although there are cusps at both B and C, large amplitudes are 
to be expected at the distance corresponding to C, but not to B. 


3- Second model 


In the second model, we take circumstances as in the first model, except 
that there is now a discontinuous jump in v at the level 7}. Let mq and y», where 
N» < a, be the values of 7 just above and just below this level. 
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The time-curve still has the appearance of Figure 3, but the branch BC now 
corresponds to rays which lie entirely in L and are reflected upward at the boun- 
dary r = 7. 


For a < p < % (the branch AB; unreflected rays in L), (3), (4) and (5) 
continue to apply. 


For 9» < p < %a (the branch BC; reflected rays), we have 
A = ag{cos“(p/no) — cos1(p/na)}, (10) 
t = ao{(no*— p*)t — (na? —p*)}}, (11) 
dA/dp =a0{ — (no? — p?) + + (na? —p?) +}. (12) 




















FIG. 3. 


For p < 7 (the branch CD; rays entering N), 


A = ap{cos~1(p/n0) — cos(p/na) + « cos(p/nv)}, (13) 

t = ao{(no?—p*)! — (na — p)* + «(m0? — p*)*}, (14) 

dA|dp = ao{—(n0*—p*)* + (na? —p?)* — (qo? — p?) 4}. (15) 

These equations show that A and ¢ are again continuous throughout the whole 


range of values of p, but that dA/dp is discontinuous, and in general not zero, at 
both B and C. 
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Hence in this case, large wave amplitudes are to be expected at neither of the 
distances corresponding to B and C. 


4- Third model 


In the first two models, there were discontinuous changes in 9 or « between 
L and N. We now consider a third model, containing L as before, but between 
L and N containing a layer M, for which re < r < 7 and in which « decreases 
continuously from ap to kao. (We continue to take « constant and equal to a, 
xag in L, N, respectively, with x < 1; and we take 7 continuous throughout L, 
M and N; N now occupies the region r < re.) 

For mathematical convenience, we take 


da/dn = cn (16) 
in M, where c is a positive constant, so that 
2a(1— x) = ¢(m?—72"). (17) 


It can be deduced (from equations in the writer’s text (1959, Section 7.23)) 
that, in general, 


dA ag 20 du/dy 
— = —- —— + Mn. 
dp V(m0°—P*) 3s VC? -P*) 





(18) 


For rays entirely in L, (18) yields (5), as usual. For rays with lowest points in M, 
(18) yields 


satan ~ > + 
dp Vf") 





dA ng { *  da/dy 
Vv 


. an 
(4? —p*) 


— axolno?— p2)-+ + e(m?—p®), (*9) 


by (16). And for rays which enter N, 


dA|dp = —aolmo?— p*)-+ + e(m®—p*)'— o(ma*— p*). (20) 


By (5), (19) and (20), dA/dp is now continuous for the whole range of values 
of p concerned. For rays entirely in L, dA/dp is negative (as usual). For rays enter- 
ing M, it can be deduced from (19) that, as p decreases, dA/dp passes continuously 
through a single zero (corresponding to B” in Figure 4) in the range y2 < p < 7, 
and so becomes positive, provided 


¢ > ao{(q? — n2”)(no? — n2?)} +, 
1.e., by (16), provided 
K < 1—4(m?—2?)*(n0? — 27). (21) 


For rays penetrating N, as p further decreases, dA/dp may change back from 
positive to negative, passing through a second zero (the point C” of Figure 4). 
This can be inferred from the fact that, for p = 0, (20) gives 


d\/dp = —aono*+¢(m — 72); (22) 
and, by (16), the right side of (22) is negative if 


K> 1—(1+72)/20. (23) 
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There is no difficulty in finding values of 7/n9 and 2/70 which satisfy both the 
conditions (21) and (23) with « < 1, and which give sensible models. Thus the 
curve shown in Figure 4 gives a further type of possible behaviour of dA/dp. 
The corresponding p—A curve is shown in Figure 5, while the t—A curve has the 
same general appearance as Figure 3, with triplication again. 

The interesting feature of the present model is that dA/dp is now zero for both 
B and C, so that large amplitudes may occur at the distances corresponding to 
both cusps. 


5. Discussion 


The models presented show that, when there is triplication in the time-curve 
as in Figure 3, three main cases arise, in which there may be abnormal wave 
amplitudes at distances corresponding (a) to neither B nor C, (6) to C but not B, 
(c) to both B and C. 














Fic. 4. Fic. 5. 


The discrimination between (a) and (b) rests on whether v is continuous or 
not. Even a small jump in v at the appropriate level can substantially reduce 
abnormally large amplitudes that might otherwise occur. 

The discrimination between (b) and (c) rests on the fine distinction between a 
sudden and a rapid but continuous change in the gradient —dv/dr. 

The analysis shows that sudden changes in the second or higher derivatives of 
v with respect to r, as well as sudden changes in the first derivative or in v itself, 
can have important influences on amplitudes. Thus further illustration is pro- 
vided of the danger in trying to infer even broad features of velocity variations 
from amplitude data. 

A central result is that, in the important case of a so-called ‘‘second-order” 
discontinuity in v, i.e. a sudden increase in —dv/dr or decrease in «, large ampli- 
tudes are not in general associated with the cusp B. On the other hand, large 
amplitudes may be associated with B if the discontinuity is what might be called a 
“third-order” one. 

The analysis for the case in which wv suffers a sudden diminution with increase 
of depth is being set down elsewhere; this case also lends itself readily to dis- 
cussion with the use of 7 and a. It is easy to show that the cusp which arises 
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in this case (see e.g. Bullen 1959, p. 117, Figure 28) is in general associated with 
large amplitudes, corresponding for example to the well-known observations of 
PKP waves near 142°. 
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Note on the Paper of A. H. Cook, ‘“‘The External Gravity 
Field of a Rotating Spheroid to the order of e*’’ 


Walter D. Lambert 


(Received 1960 March 14) 


Summary 


A. H. Cook published in the September 1959 number of this Journal 
an article entitled ‘““The External Gravity Field of a Rotating Spheroid 
to the order of e8”’. Cook’s formulas for the potential of such a spheroid 
stop with the term in P¢(sin¢), which is given, as indicated, by the 
title, to an accuracy of e°. This no’e states, without proof, formulas 
for the coefficients of P2,(sin ¢) with no limit on nm and no limit to the 
formal numerical accuracy of such coefficients. 


Some of Cook’s other formulas may be simplified by a slight change - 
in notation. 


Dr Cook’s paper appeared in the September 1959 number of this Journal, 
pp. 199-214. It should be noted that the spheroid in question is an exact ellipsoid 
of revolution, assumed to be an exact surface of equilibrium under the combined 
effects of its own gravitational attraction and the centrifugal effects of uniform 
rotation about its minor axis. 

Cook gives, formula (5.2), page 211, formulas for the J’s, the coefficients of the 
expansion of the gravitational potential of the ellipsoid in terms of the zonal har- 
monics of the geocentric latitude multiplied by appropriate powers of the reciprocal 
of the radius vector. Cook’s formulas are accurate to the third order of small 
quantities, such as the flattening, and include the zonal harmonic of degree six. 
The purpose of this note is to show that perfectly general formulas, accurate to 
any desired order of smallness, may be given for terms of any degree. Also, some 
of Cook’s other formulas can be slightly simplified by using a slightly different 
parameter to express the centrifugal effect. Only the results are stated here. The 
proofs are given in Technical Paper (716)-2-261 of the Mapping and Charting 
Research Laboratory of The Ohio State University Research Foundation, Colum- 
bus, Ohio. 

For various reasons, it seems inconvenient to adhere entirely to Cook’s nota- 
tion. Instead, the following symbols will be used: 
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Semi-major axis of ellipsoid a a 
Semi-minor axis of ellipsoid b b 
First eccentricity of ellipsoid 

Second eccentricity of ellipsoid 


Flattening of ellipsoid - e 


Mass of ellipsoid M 
Newtonian constant of gravitation f 
Radius vector of point in space r 
Geocentric latitude of point in space 6 
Angular velocity of rotation of ellipsoid w 
Zonal harmonic of degree 2n P2p(sin 6) P2»p(sin 6) 


. a\2n ‘ , 
Coefficient of | *) P2,(sin @) in the 
r 


expansion of the potential Aon Jon 
Mass of the spheroid M M 
Other symbols will be introduced as needed. 


The potential of the ellipsoid under the assumed conditions can be expressed 
in closed form in terms of curvilinear coordinates. This has been done more than 
once. But curvilinear coordinates are much less convenient in practice than the 
ordinary spherical coordinates (r,@). Because of symmetry no third coordinate 
(corresponding to longitude) is needed. ‘The direct transformation of the formulas 
for curvilinear coordinates into the corresponding formulas for spherical co- 
ordinates is laborious and leads to no general results. This difficulty can be 
avoided, as here, by the use of a highly artificial scheme of densities: A homo- 
geneous fluid ellipsoid in equilibrium with a superficial coating of high negative 
density. We know that a given gravitational field can be produced by infinitely 
many distributions of density and, since this artificial scheme of densities can be 
shown to give the same formula for gravity on the surface of the ellipsoid as the 
formula derived from curvilinear coordinates, which is independent of any special 
scheme of density, the results for the artificial density scheme are valid also for 
exterior space. 

The formulas in terms of curvilinear coordinates involve the use of Legendre 
functions of the second kind with an imaginary variable icot x«. These are needed, 
incidentally, in the present discussion. We assume 


¥ tan-le 
gia) = —Q,(tcota) = 1—acota = 1— » 
€ 


goa) = —iQde(icot a) : (3 cot?« + 1)x—3 cot a| 
2L 


2 
(° si ) an-e—3]. 
| e2 € 





We have also 


dqo(x) 


= 3 .cosec*ag;(«)—1. 
da 
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For small values of « or « these forms are inconvenient for practical computa- 
tion, and expressions as power series are preferable. The series are often expressible 
as hypergeometric series, a fact that makes it easy to pass from a series in € to a 
series in e. In the customary notation 


g2(2) = —3F (5 im -) 
wz \s 2 


2 . .. - 
= <al1 — —€2 + 4 —-—<6 + | 
15 + tt ® 


= < el: $8 ty Dees | -|. 


15 14 56 528 


In what follows it is convenient to use a modification of g2, namely, 4 = 3q2/e, 
or in series form 


2 
——e4 + —<6 ——-<8 + ... 


~~ Ft BB 


I 155 
2 Big 38 oy 
35 231 


We shall also need the quantity 4; and its reciprocal; 


n= 9/(=). 


8 
14-240 . A—-—... 
7 231 


I 64 


1+ a ee mA 
7 231 
cy Ven Wie 
7 49 11319 
8 6 
tofe— figs... ©) 
7 49 11319 


Imagine an unyielding sphere of radius a and uniform density p;. Gravity on 
its surface will be 


= rhe. (7) 


Let this sphere rotate with the same angular velocity as the Earth. The angular 


acceleration at the equator of rotation will be w2a. Let us form the ratio of this 
to gravity G and call it m: 


wa Zw 


my) = . 
4mkpy 
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Now instead of a sphere of unyielding material let us suppose that we have a 
fluid ellipsoid of the same density p; throughout and with an equatorial radius 
equal to a. Its mass, because of the flattening, will be slightly less than that of 
the sphere of radius a. The condition for equilibrium of this fluid ellipsoid is 


(8) 


If p; were the density of our actual Earth, this quantity m would be equal to the 
quantity that Cook calls m. See middle of page 204 after equation (3.8). 

But if we take p; as the density of the actual Earth, the flattening comes out 
approximately equal to jm. This is too large. Instead let us take the flattening 
what we know it to be, say 1/297 for the International Ellipsoid. ‘This means that 
pi comes out as greater than 7 g/cm®, which exceeds the mean density of the Earth. 
To bring the mean density down, let us apply a surface coating of sufficiently high 
negative density, a coating such as was used long ago by Pizzetti. However 
absurd this set-up may be in a physical sense, the mathematical conclusions derived 
from its use are sound enough. 

The potential V; of the homogeneous ellipsoid of density p;, and eccentricity e 
for an exterior point whose coordinates are (r, @) is given by 


47ka®b (my I 1 /ae\? 1 /ae\4 
nea BYLE (Syne HE) 
r $/Lr3 3°5\ 7. 57\ 7 


\° (-1" 
ry, ovkiiaaiad (2n+ ee FH Pont |. 


The argument sin@ is understood with the P’s, which represent, as usual, zonal 
harmonics. 

We now use the fact that the potential of a homoeoid is uniform over the 
surface. Let pda be the surface density of the coating, whose mass is to be sub- 
tracted from that of the too massive fluid ellipsoid to bring the total mass down 
to that of the Earth. We understand p to be very large and da to be very small, 
with their product approaching the necessary finite limit, the total mass being 


(9) 





M= <li, (10) 
3 


where pm is the mean density of the Earth. The limiting value of pda turns out 
to be 


pda a (= 

= —Ppm|——1}. 
3 \ yp 

This gives for V2, the potential of the coating 


kab \2 4 
V4 = S(t) net) ny 
3 °f ys 3\r s\r Nr 


+ ai? is (=) "Pea + mt 


an+iI\ r 
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Substitute for % its equal Ze%,. The resultant potential V = Vi—V2 may be 
written 


/ ‘a\2 I 10 my, \ /a\* 
Fah~me Pet iet(e-— (*) P.- 
\ /\p, 5 7 vn 
~*a(2_37)/(* 
7. $a 


(—s) 
cla eae ) "Peet 


By using expression (6) for 1/1 it is seen that equation (13) checks perfectly 
with Cook’s equation (5.2), page 211, where the various coefficients Jon (here 
called Agn) are expressed to quantities of the third order of smallness. Formula 
(13), however, enables any coefficient, even those beyond Ag, to be computed to 
any desired degrees of formal accuracy. For the present, however, this advantage 
is theoretical and esthetic rather than practical. Cook remarks, in effect, that it 
seems more logical and more convenient to take the coefficient (a/r)?"Pon as it 
stands rather than to tie it up to a numerical coefficient, the reciprocal of the 
highest power of sin @ in P,(sin@). This is what is done in the JD notation, where 


2 8 
A2 = —J and Ay = —D. 
3 35 


A similar opinion as to the most convenient notation has been expressed by 
Brouwer (1959). 


The quantity 


Yaate) 


is positive up to m = 5, but for m = 5 it is extremely small and for m = 6 and 
beyond it is negative. 

The advantage of using ys; rather than ¢ is that 4 is very nearly equal to unity, 
so that my, (m/:1) and m (to be defined later) differ from one another by quantities 
of the second order of smallness and may be equated to one another in approxi- 
mate comparisons of different expressions for the same quantity. 

Cook’s m, which is the m here used, differs slightly from what may be called 
the traditional m, the ratio of the centrifugal acceleration at the equator to gravity 
at the equator. Denote equatorial gravity by ge. Then 


(14) 


For an exact relation between m and mj, as here used, we define a new function 7 
based on expressions (1) by the equations 
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or in terms of the hypergeometric function 


n= 5@F(1,2,4-2 
Sa ae fleas 


3 2 8 16 
‘ 5 eae sane 
5 7 63 231 


For the present we need the quotient »/) rather than 7 alone 


By q’2 is meant the derivative of g2 with respect to «. By division weffind 


F 6 208 
. 1+n8- t+ a, 
op L 87 «6441 11319 


j I oll 
: r+at4 ay a. 


, oo 33 957 





r 3. = 068 
: 1+—f+—2fe+- f...| 
2 


L 7 44! 3773 
The relation between m and m, as here defined, is 


m 
my) 


] 
I+m- 


my 


7 
[a= 


To small quantities of the third order we have 


3 
m, = wn — on? — Fant on, (18) 
2 7 


Apply this to Cook’s formula (3.9), page 204, changing the m’s to my’s to conform 
to the notation of this paper. The coefficient of sin?¢ in terms of m, and f is 


I I I 
Sy —f+ “Sony? wot m — —f*m, fin + m3. 
2 4 14 4 8 


Replace m, by m. There results 
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Three terms have cancelled out and the remaining third-order term has an un- 
usually small coefficient. The coefficient of sin?2¢ is changed from 


I 5 17 15 
ee ee! + aes 
gf gmt son gms 


<f 2 =mf + nt (20) 


a cancellation of one third-order term. The coefficient of sin*¢sin?2¢, being 
already of the third order, is not simplified by the change. 


The quantity / can be made to appear in Cook’s equation (3.8a), page 204. 
That equation is equivalent to 


— = ge 


kM 
ab” &| 


+m), (21) 


ip 
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The Palaeomagnetic Poles for the Lower Jurassic of 
Europe 


K. M. Creer,* E. Irvingt and S. K. Runcorn* 


(Received 1960 March 31) 


From his palaeomagnetic study of some Lower Jurassic rocks of N.W. Europe, 
Girdler (1959) concluded that the direction of the geomagnetic field was signifi- 
cantly different from the present dipole field, corresponding to a pole lying in the 
region of the Caspian Basin, a result which is quite different from that inferred 
from the polar wandering diagram for Europe given by Creer, Irving and Runcorn 
(1957) and different from that obtained by Hargraves and Fischer (1959) from 
some Lower and Middle Jurassic beds of Central Europe. But we consider that 
there is an alternative interpretation of his results. 

The directions of magnetization of many of his formations, the Bridport sands, 
Yeovil sands, normally magnetized specimens of the Cotswold sands, yellow sands 
and grey sands of the Blea Wyke Beds, stereograms of which are shown in Figures 
2, 3, 5b, 6a and 6b respectively of his paper, are grouped around the present dipole 
field—a result which shows either that the geomagnetic field in Lower Jurassic 
times was similar in direction to the present dipole field, or that, if it was not, the 
original magnetization decayed and was replaced by one acquired in recent times. 
Girdler presents no evidence but on page 355 of his paper he presumes the latter. 
His conclusion, therefore, so far as his English work is concerned is based only on 
the directions of magnetization of the Midford sands and the reversely magnetized 
specimens of the Cotswold sands (Figs. 4a and §a respectively of his paper), which 
being different from the present dipole field are thought by him to be primary. In 
fact these directions are scattered sufficiently close to a plane passing through 
the present dipole field, as is shown in the polar equal area projections below, to 
cause one to suspect the presence of two magnetizations, (a) the primary or original 
one and (b) a secondary or later one along the present dipole field. Girdler 
does not mention this possibility, although it has been demonstrated in sediments 
many times before, e.g. by Kawai and Kume (1953), by Creer (1955, 1957) and 
Runcorn (1956). Such a secondary component is stable in the sense that repeat 
measurements after six months or more show no change in its direction or inten- 
sity. A well-known test for stability (Graham’s (1949) fold test) is applied by 
Girdler to the Midford sands, in which he says that the directions of magnetiza- 
tion are less scattered after correction is made to the geological dip of the beds 
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than before (see Figs. 4b and 4a respectively of his paper), although streaking is 
still present. 

Figure 1 combines the data on Girdler’s Figures 4a and 4b; the dots are direc- 
tions of magnetization of specimens with respect to the present horizontal, but 
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those from tilted beds are joined by lines to the directions of magnetization, repre- 
sented by crosses, which are corrected for geological dip. These two latter sets of 
directions lie along a great circle passing through the present dipole field, repre- 
sented by a full square. Those directions in the N.W. quadrant of the diagram, 
which have not been corrected for geological dip, are likely to be the vector resul- 
tant directions of asecondary magnetizationalong the present dipole field superposed 
on a magnetization along the Jurassic geomagnetic field, whichis not very different in 
direction from the present one. This magnetization was apparently acquired before 
the rocks were tilted to their present dip, a geological tilt which would cause this 
primary magnetization to move to a position with respect to present horizontal, 
indicated by the large open circle on the figure. Thus we do not consider that 
Girdler’s argument for the stability of these rocks carries any weight and, of course, 
no arguments for the stability of the directions of magnetization in flat lying beds 
in the Midford sands or the “reversely magnetized” specimens of the Cotswold 
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sands shown in Figure 2 are available. They of course may, too, be suspected of 
secondary magnetization, but if so, the primary magnetization must be assumed to 
be one opposed to the present dipole field, represented by the open square in 


N 


» 
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Figures 1 and 2. It is however clear that these points are not well fitted by any 
one great circle, but as there are an infinity of such joining the supposed primary 
and secondary magnetizations of these rocks, this may not be surprising. Alter- 
natively, the divergences from this plane may be due, as Creer (1957) has shown, 
to a temporary component of magnetization present in the specimens acquired 
during the experiments. 

Girdler supports his conclusion concerning the direction of the geomagnetic 
field in Jurassic times in England by measurements on some tuffs of Hettangian 
age from the northern Pyrenees. Again Graham’s fold test is applied and it is 
stated, although it is by no means obvious, that the directions of magnetization are 
less scattered after corrections have been made for geological dip (Figure 8b) than 
before (Figure 8a). The use of a statistical method appropriate to a symmetrical 
distribution is hardly appropriate in this case. Further, Girdler’s statement that 
“when geological corrections are applied the directions of magnetization become 
less scattered” is even more surprising as he states that through lack of adequate 
field mapping the assumed geological dip is subject to an error of + 15° (p. 360). 
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Finally, it may be noted that the directions of magnetization with respect to the 
present horizontal are distributed close to a plane including the present dipole 


field (see Figure 3); this again implies the presence of a secondary component of 
magnetization. 
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\ 
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We consider that Girdler’s observations do not enable any significant state- 


ment to be made concerning the position of the geomagnetic pole relative to 
Western Europe during the Lower Jurassic. 
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The palaeomagnetic poles for the Lower Jurassic of Europe 
Comment by R. W. Girdler. 


The author is aware of a possible secondary component of magnetization 
being present in the Upper Lias Sands; this is discussed on page 355 of his paper 
and Figure 7 on page 358 was drawn to emphasize the varying degrees of stability. 

The objection concerning the Midford Sands arises because it so happens that 
applying the geological correction required for some of the samples involves 
moving the remanent vectors along great circles which pass close to the Earth’s 
present dipole field. The alternative explanation of the results given in the above 
letter is unlikely for the following reasons: 

(a) The samples not requiring geological correction do not lie close to a great 
circle passing through the Earth’s present dipole field as can been seen from 
Figure 1 of the above letter. The sands are remarkably similar (see the geological 
references cited at the end of the paper) and come from a small area in Somerset; 
it would be surprising if only those requiring geological correction should lie on 
the proposed great circle when the majority of the others do not. 

(+) The following table shows the results of calculating the mean directions of 
magnetization of the specimens requiring geological correction separately : 


6 d K N 
(a) Samples not requiring correction 97°5 69:1 . 8-8 30 
(b) Geologically corrected samples 118°1 70°6 11'S 15°3 12 


(c) Mean of all samples 103'5 698 7°4 ° 42 


(The symbols have the same meaning as for the table given on page 358 of the 
author’s paper.) 

‘The mean directions together with their circles of confidence at the 95 per cent 
level are illustrated in Figure 1. These statistics are calculated from the geologi- 
cally corrected results which do not show obvious streaking along great circles 
through the present dipole field. The agreement for specimens requiring geological 
corrections and those not requiring geological corrections is very good and would 
hardly be expected if a secondary component of magnetization is important. 

(c) The mean direction of all the samples is approximately 180° different from 
that of the reversely magnetized Cotswold Sands which are of nearly the same age 
(Girdler 1959). Previously (e.g. Blackett (1956), Irving & Green (1957) and others), 
this has been considered evidence of stability. It is probably significant that the 
writers of the letter note that the reversely magnetized Cotswold Sands ‘‘are not 
well fitted by any one great circle’’. 

It is concluded that there is no good evidence, unlike the excellent example 
of Creer (1957), that a secondary component of magnetization parallel to the 
present dipole field is important in the Midford and Cotswold Sands. 

Throughout this work much thought has been given to stability and it is unfor- 
tunate that the sediments are too weakly magnetized to enable laboratory experi- 
ments to be performed. The Upper Lias Sands of the west of England were 
deposited continuously as one unit in a time interval of probably less than 
6 x 108 years and this is comparable, if not less than the time covered in sampling 
by other workers for other periods. As previously pointed out, some of the Upper 
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Lias Sands do show a secondary component of magnetization along the direction 
of the present dipole field. Figure 7 has been redrawn (Figure 2) with the reversely 
magnetized Cotswold Sands plotted as a normal direction. The great circle 
through the present dipole field and the mean directions of magnetization is shown 
as a dotted line. In the absence of any other data the directions of n.r.m. of the 
Cotswold and Midford Sands (180° opposed) are taken to be close to the field 
direction for the Jurassic period. 


Fic. 1.—Mean directions of natural remanent magnetization of Midford 
and Cotswold Sands together with circles of confidence at the 95 per cent 
level. (a) For samples not requiring geological correction, (6) for geo- 
logically corrected samples, and (c) for all samples. The mean direction 
for the Cotswold Sands (d) is shown with the circle of confidence as a 


dashed line. 


‘The work was an attempt to fill a gap in our knowledge of palaeomagnetism, 
very few results having been obtained for the Jurassic period. Most interesting 
is the presence of an easterly declination (slightly larger than for the Triassic) 
and the fairly steep inclination of 67°. The writers of the letter make no comment 
on the fact that the latter is very close to that found by Nairn (1957), discussed 
on p. 359, and also on the fact that the rocks of Hargraves & Fischer come from 
a thrust area and are unusually tightly clustered. Finally, it is perhaps not sur- 
prising that the result does not agree with the smoothed polar wander curve of 
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Creer, Irving & Runcorn for this assumes no continental movements or local 
movements of the Earth’s crust and as more and more results become available 
it is possible that their polar wander curve will become much more complex. 


R. W. GIRDLER 
Durham Colleges, 
Science Laboratories, 
South Road, 
Durham City: 
1960 June 20, 
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Fic. 2.—Mean directions of natural remanent magnetization of all the 
Upper Lias Sands with great circle passing through the present dipole 
field. 
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Geophysical Discussion of the Royal Astronomical 
Society 


A discussion of papers received during the session was held in the rooms of the 
Royal Astronomical Society on 1960 April 22 at 16h 15m under the Chairmanship 
of Dr R. Stoneley, F.R.S. 

Dr K. F. Bowden, of the Tidal Institute, Liverpool, presented a paper written 
by himself, L. A. Fairbairn and P. Hughes, on The Distribution of Shearing 
Stresses in a Tidal Current (1959). The shearing stress in a uniform liquid 
flowing steadily along an open channel, increased linearly from the surface to the 
bottom, regardless of turbulence, and the distribution of velocity with depth could 
be determined theoretically. When the flow varied with time, and effects due to 
the Earth’s rotation were added, as in the case of tidal flow, no direct inference 
could be drawn about the distribution of shearing stress with depth. ‘The investi- 
gation discussed was aimed at determining the shearing stress as a function of 
depth from measurements of the vertical distribution of current throughout the 
tidal cycle and from estimates of the bottom friction. The equations of the system 
were given. Current measurements were made in Red Wharf Bay, Anglesey, 
North Wales, from an anchored vessel, using a Doodson recording current meter 
augmented by a Charnock meter for some near-bottom measurements. ‘The velo- 
city profiles were determined and the bottom stress was calculated. The stress 
determination gave a nearly linear profile from surface to bottom at the times of 
maximum flood and maximum ebb, but departed strongly at other stages of the 
tidal cycle. The eddy viscosity tended to be higher at mid-depths than at the 
boundaries and to reach a maximum when the current was greatest, about three 
hours before and after high water. 

In the discussion, the Chairman asked what was the effect on the stress deter- 
minations of viscous drag due to wind on the sea’s surface. The speaker answered 
that, with a 10-knot wind, the stress on the surface was o-6dyn/cm? compared 
with a figure between 3 and 4dyn/cm? at the bottom and had been neglected. 
It would have been an improvement to take this into account but would have 
required measuring and recording the wind speeds, and they had not the instru- 
ments for this. Professor Sheppard asked if the physical causes were clear as to 
why the eddy viscosity was different between the flood and ebb of the tide. 
Dr Bowden said the case was not clearly understood; there was no density 
gradient to account for it, but the current velocity was slightly different in 
the two cases. 

Dr K. M. Creer, of King’s College, Newcastle, followed and delivered the 
paper written by himself, E. Irving and A. Nairn on Palaeomagnetism of the 
Great Whin Sill (1959). This described the results of a study of the remanent 
magnetism in rock samples collected at thirty-two sites in the Whin Sill of 
Northumberland and Durham. The authors also reported that the mean direc- 
tions of magnetization had an oval distribution, discussed the possibility that this 
scatter reflected the geomagnetic secular variation during the time of interjection, 
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and proposed a method for calculating the errors associated with the palaeo- 
magnetic pole positions. It was concluded that a direction of 187°-8 declination 
and —4°-9 inclination represented the vector mean of the magnetizations they 
had observed. ‘There was not time for discussion of this paper. 

The third paper was by Mr J. E. Jackson, of Cambridge University, on 
Observations with the Cambridge Pendulum Apparatus in North, Central 
and South America in 1958 (1959). In introducing this report the author 
reviewed the needs for geodetic study which prompted these observations. About 
100000 stations between which the gravity differences should be known to a 
1 mgal accuracy were required for accurate knowledge of the shape of the Earth, 
these should be fairly uniformly distributed over the Earth's surface, and at 
several stations the absolute value of gravity should be measured. ‘The Cambridge 
Pendulum Apparatus was used in making the observations reported. ‘This con- 
sisted of four invar pendulums swung in pairs. Eleven stations were occupied 
both on the outward and on the return journeys from England to South America: 
‘Teddington, Madison, Mexico City, Panama, Caracas, Quito, Lima, La Paz, 
Santiago, Buenos Aires and Rio de Janeiro. ‘The computed gravity differences 
were given. Several slides were shown of the gravity stations and their surround- 
ings. In the discussion which followed Dr A. H. Cook remarked there were now 
five stations where absolute measurements of g had been made. 

Dr J. C. Weston of the Department of Scientific and Industrial Research spoke 
on his paper, Underground Explosions (1960). ‘The low frequency of the energy 
radiated from underground explosions was proportional to the square of the fre- 
quency of the seismic wave (f?) and to the square of the total energy (W2). The 
f* law was derived from relationships between pulse shapes and the radiated 
spectrum. ‘The W? law came from a consideration of body waves. Underwater 
explosions had provided evidence for these laws, but in underground explosions 
only the W2 law had been verified. Some comments were made about the fre- 
quency distribution in underground and underwater explosions, and about the 
theoretical source level and critical radiation in underground explosions. ‘There 
was no discussion. 

Sir Harold Jeffreys was the last speaker and delivered a paper by Professor 
K. E. Bullen of Sydney University on The Derivation of Seismic Velocities 
at Depth from Travel-Time Curves (1960). In deriving the distribution of 
velocity with depth from travel-times, the usual method is to find dt/dA, and then 
an integral depending on dt/dA had to be worked out numerically to give a relation 
between A and the maximum depth reached by the ray. In practice the work was 
rather laborious. What made it worse was that one might want to examine the 
effect of changing f over some range of A, and the only way of finding it was to do 
all the work again. It was known that there was an exact solution of w of the form 
a— fr, but Bullen pointed out some time ago that there was another, if v were of 
the form Ar-*, when the times were proportional to sing(k+1)A. The paper 
depended on using this as a first approximation. With Bullen’s treatment a first 
approximation did facilitate later steps. He applied it to the Earth’s main core 
and got results agreeing closely with those found by direct integration. Jeffreys 
commented that it was rather odd that it should be so because d2c/dr? is positive 
for the first approximation and the actual values were usually negative, as for 
v = «a—Br?, but the method did manage it. Effectively it introduced d loge/d logr 
and when this was found c came by one integration of a fairly small quantity. 

The Chairman thanked the speakers and adjourned the discussion at 18h 10m. 





376 Geophysical discussion of the Royal Astronomical Society 


References 


Bowden, K. F., Fairbairn, L. A. & Hughes, D., 1959. Geophys. J., 2, 288. 
Bullen, K. E., 1960. Geophys. J., 3, 258. 

Creer, K. M., Irving, E. & Nairn, A. E. M., 1959. Geophys. J., 2, 306. 
Jackson, J. E., 1959. Geophys. J., 2, 337. 

Weston, D. E., 1960. Geophys. J., 3, 191. 





Book Reviews 


Palaeomagnetism and Stratigraphic Correlation 


A. N. Khramov 


(published by the Government Science publishers of petroleum and mineral 
fuel literature, Leningrad 1958, 204 pp. In Russian.) 


This book, its author says, was written “because of the unbearable absence in 
Russia of a book explaining the present state of knowledge and methods of palaeo- 
magnetism”’, and “‘is directed at a wide circle of geologists and geophysicists and 
also at the specialist who is beginning palaeomagnetic research with an eye to its 
use both in stratigraphical problems and in tracing the geological history of the 
Earth’s magnetic pole’. 

It is divided into three chapters. The first contains a general account of palaeo- 
magnetism and its implications for geology and geophysics based chiefly on work 
prior to 1955 in countries other than Russia. ‘There are many small inaccuracies 
in reporting and also some curious referencing, for instance, Table 7, page 48, 
refers to Hospers (1955) for studies of the lava flows of Mt Etna which were made 
in 1926 by Chavallier, but it is likely that Mr Khramov has preferred to give 
reference to reviews in journals available in Russian translation when the originals 
are not available. The second chapter deals with palaeomagnetic work in ‘Turk- 
menia where beds from Jurassic to Pleistocene have been studied, and the third 
deals with work on Ordovician to Lower Triassic rocks from the Russian Plat- 
form. The second chapter should interest oil geologists and the second and third 
chapters are necessary reading for research workers in palaeomagnetism. 

The Russian work confirms the two most important features of palaeomagnetic 
results from other countries, namely the occurrence of serial reversals of polarity 
superimposed on the steady drift of the directions from period to period. Mr 
Khramov, following most ‘“‘western” workers, interprets the former as due to 
reversals of the main geomagnetic field and the latter as due to polar wandering. 
The discrepancy between the pole results obtained from different continents he 
regards as strong evidence of continental drift and bemoans the unnecessary 
hostility to this idea amongst Soviet geologists. 

The work in Turkmenia arose from the need for more accurate stratigraphic 
correlation of rock formations whose equivalents in Caucasia yield oil. A close 
study of Upper Tertiary and Quaternary beds has been made with less detailed 
work on Jurassic, Cretaceous and Lower Tertiary formations. An alternating 
succession of “normal” and “reversed” zones has been established from surface 
sections, and this is correlated with similar effects in bore-core material. Several 
promising methods of distinguishing between different zones of the same polarity 
are proposed and applied with some success. Results from transition beds between 
“normal” and “reversed” zones suggest that reversal of the Earth’s field is not 
due to decay of the Earth’s dipole followed by build up in the opposite direction, 
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but rather to its toppling over through 180°. Regarding the pole positions obtained 
from the Turkmenian work it is of interest that the north pole for the Plio-Pleisto- 
cene has the co-ordinates 75 N, 122 W. ‘This revives the opinion once held that a 
polar displacement towards North America was necessary to explain the assymetry 
of the Pleistocene glaciation. Similar effects should, but have not turned up in 
work elsewhere, but this could be due to inadequate correlations and the limited 
data. 

The work in European Russia is notable chiefly for studies of the Permian 
classical sections in the eastern Urals. Except for the highest levels of the Permian 
(the ‘Tartarian which could in part be Lower Triassic) the directions are all 
“reversed”’ as for Permian beds elsewhere in the world. ‘The pole agrees very well 
with that obtained from the western European Permian. The level of the return 
to “normal” polarity in the Tartarian promises to be a useful world-wide datum. 

A translation into English of this book will be published in 1960 or early 1961 
by the Australian National University. 


Introduction to Geophysical Prospecting 
M. B. Dobrin 
(2nd edition, 1960, 446 pp., McGraw-Hill, 74s.) 


This is the second edition of the book originally published in 1952. It is a 
second edition in more than name: new sections have been added at the expense 
of some of those in the first edition and the text with its very useful references has 
been thoroughly revised to bring it up to date. The number of revisions illus- 
trates the appreciable extent to which prospecting techniques have been modified 
in the past decade. ‘The basic methods are unchanged but there have been con- 
siderable developments in both the instruments and the interpretation techniques. 
Particularly in seismic work, these have coincidentally led to a better appreciation 
of the effect of the recording process on the data quality. 

Seismic methods occupy a predominant position in the industry at the present 
day and the sequence of chapters in this edition has been revised to emphasize 
this, with additional space being devoted to both reflection and refraction tech- 
niques. ‘he new material here includes magnetic tape recording, corrected record 
sections, geophone patterns and recent interpretation methods. One particu- 
larly interesting addition is a chapter on the geological interpretation of reflection 
records with its discussion of geological/geophysical coordination. ‘The author is 
to be congratulated on avoiding the woolly talk which usually surrounds this topic 
and on illustrating in practical terms why this coordination is required. 

Developments in gravity and magnetic surveys have been less spectacular than 
those in seismic surveys but the sections in interpretation have been extended and 
the latest instrumental developments included. The chapters on radioactive and 
electrical methods have also been enlarged although in the case of electrical methods 
perhaps not sufficiently to satisfy the mining geophysicist. This extra material 
has been incorporated at the expense of three chapters of the first edition: of these, 
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the only serious loss is that dealing with well-logging methods although its 
omission is justifiable in view of the tendency of well-logging to become an 
independent geophysical field. 

The short chapter on earthquake seismology and the structure of the Earth 
has been retained; in spite of the author’s apologia for its presence, it seems some- 
what out of place in this book. Whilst commercial geophysics owes a great debt 
to earthquake seismology, one would hope that any geophysics student could be 
left to acquire this background information elsewhere. The least satisfactory part 
of the book is perhaps that dealing with seismic instruments. Although it includes 
recent developments and has been largely rewritten, there are a number of minor 
errors and rather misleading statements. 

In spite of these minor criticisms, this book, like its first edition, is very com- 
prehensive, well-produced and very readable. It can be firmly recommended 
to its intended readership of students and of those practising geophysicists who 
require a quick reference to fields outside their own speciality. 


A.T.D. 
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